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■ We review the status of covariant methods in quantum field theory and 

' quantum gravity, in particular, some recent progress in the calculation of 

the effective action via the heat kernel method. We study the heat ker- 
nel associated with an elliptic second-order partial differential operator of 
Laplace type acting on smooth sections of a vector bundle over a Rie- 
mannian manifold without boundary. We develop a manifestly covariant 
method for computation of the heat kernel asymptotic expansion as well 
as new algebraic methods for calculation of the heat kernel for covariantly 
constant background, in particular, on homogeneous bundles over symmet- 
ric spaces, which enables one to compute the low-energy non-perturbative 
effective action. 
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1 Introduction 

One of the most important problems of modem fundamental physics is the prob- 
lem of reconciling classical general relativity, the theory of macroscopic gravi- 
tational phenomena, with quantum theory, so-called Quantum Gravity problem. 
This is a really difficult task since one has to answer the very basic questions con- 
cerning the local and the global structure of the spacetime itself as well as deep 
questions about the nature of quantum mechanics. 

Although, over the last several decades many competing approaches (Eu- 
clidean path integrals, string theory, loop gravity, non-commutative geometry, 
asymptotic safety, various lattice approaches and others) has been put forward and 
despite some real progress in some of these approaches in the last two decades, we 
still do not have a complete consistent theory of quantum gravitational phenom- 
ena. It looks like we are missing an important piece of the puzzle which prevents 
us to find the solution. 

In this situation it seems to be wise to go back and to recall some pioneering 
works in quantum gravity. This review will concentrate on so-called covariant 
methods in quantum gravity. Some other approaches are reviewed by other lec- 
turers of this school. The basis of the covariant methods in quantum gravity is the 
background field method. This method was developed mainly by De Witt in his 
classical papers [|2D1 and reviews (221 1221 (for the latest update see the book 
Il24l ). It is a generalization of the method of generating functionals in quantum 
field theory developed and successfully used by Schwinger ||29ll30ll . For a detailed 
review see, for example, lfT8l[T6ll26ll ). 

The basic object in the background field method is the effective action. The 
effective action is a functional of the background fields that encodes, in principle, 
all the information of quantum field theory. It determines the full one-point propa- 
gator and the full vertex functions and, hence, the whole S -matrix. Moreover, the 
variation of the effective action gives the effective equations for the background 
fields, which makes it possible to study the back-reaction of quantum processes 
on the classical background. In particular, the low energy effective action (called 
the effective potential) is the most appropriate tool for investigating the structure 
of the physical vacuum in quantum field theory. 

The only practical method for the calculation of the effective action is the 
semi-classical perturbative expansion of the path integral in the number of loops. 
All fields are split in background classical parts and quantum perturbations prop- 
agating on this background and the classical action is expanded in quantum fields. 
Then the quadratic part determines the propagators of the quantum fields and the 
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higher-order terms reproduce the vertex functions of the perturbation theory. 

In the perturbation theory the effective action is expressed in terms of the prop- 
agators and the vertex functions. One of the most powerful methods to study the 
propagators is the proper time method (also called the heat kernel method, in par- 
ticular, by mathematicians), which was originally proposed by Fock [|25l and later 
generalized by Schwinger |l29l [30l| who also applied it to the calculation of the 
one-loop effective action in quantum electrodynamics. It was De Witt [EOIE^I^^ 
who perfected the proper time method; he reformulated it in the geometrical lan- 
guage and applied it to the case of gravitational field. 

At one-loop level, the contribution of the gravitational loop is of the same or- 
der as the contributions of matter fields . At low energies (lower than the Planck- 
ian energy, hc^ /G) the contribution of higher gravitational loops should be highly 
suppressed. Therefore, a semi-classical concept applies when the quantum matter 
fields together with the linearized perturbations of the gravitational field interact 
with the background gravitational field (and, probably, with the background mat- 
ter fields). This is what is usually called the one-loop quantum gravity. The main 
difficulty of quantum gravity is the fact that there is no consistent way to eliminate 
the ultraviolet divergences arising in perturbation theory, even at one loop level. 

The present review is devoted to the development of the covariant methods for 
calculation of the effective action in quantum field theory and quantum gravity. In 
Sect.2 we review the formal structure of quantum gauge field theory and quantum 
gravity and the construction of the effective action following ffTOlfTTII . In Sect. 3 
we describe the heat kernel method and develop the asymptotic expansion of the 
heat kernel following |l2l[3l|9l[IOl[IIl. In Sect. 4 we describe the local structure of 
the Green function following 0|. In Sect. 5 we develop a method for the calcula- 
tion of the heat kernel coefficients and describe their general structure following 
[l3l|9l[T0l[IIl. In Sect. 6 we compute the heat trace in the high-energy approxima- 
tion following nHHIISl- In Sect. 7 we describe our results for the calculation of 
the low-energy heat trace following our recent work [4l[51[6l|71[l2l[T3]. In Sect. 8 
we apply the obtained results to compute the low-energy one-loop effective action 
in quantum gravity. 
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2 Effective Action in Quantum Field Theory and 
Quantum Gravity 



In this section we briefly describe the standard formal construction of the gener- 
ating functional and the eff'ective action in gauge theories. The basic object of any 
physical theory is the spacetime M, which we will assume to be a n-dimensional 
manifold with the topological structure of a cylinder 



where / is an open interval of the real line (or the whole real line) and S is some 
(n- l)-dimensional manifold. The spacetime manifold is here assumed to be glob- 
ally hyperbolic and equipped with a (pseudo)-Riemannian metric g of signature 
(_ + ...+); thus, a foliation of spacetime exists into spacelike sections identical 
to S. Usually one also assumes the existence of a spin structure on M. A point 
X = (xf^) in the spacetime is described locally by the time and the space coordi- 
nates (x^,.. .,x"~^). We label the spacetime coordinates by Greek indices, which 
run from to (n - 1), and sum over repeated indices. 

Let us consider a vector bundle "V over the spacetime M each fiber of which 
is isomorphic to a vector space, V, on which the spin group Spin(l,n- 1), i.e. 
the covering group of Lorentz group, acts. The vector bundle can also have 
an additional structure on which a gauge group acts. The sections of the vector 
bundle 'V are called fields. The tensor fields describe the particles with inte- 
ger spin (bosons) while the spin-tensor fields describe particles with half-integer 
spin (fermions). Although the whole scheme can be developed for superfields (a 
combination of boson and fermion fields), we restrict ourselves in the present lec- 
ture to boson fields (which, without loss of generality can be considered real). A 
field (f is represented locally by a set of real-valued functions (p = {(p^{x)) , where 
A = 1, . . . ,dimy. Capital Latin indices will be used to label the local components 
of the fields. To construct invariant functionals we need to introduce an invariant 
fiber inner product and an inner product 



where <ivol(x) = dxg ' , g = |detg^vl> is the natural Riemannian volume element 
defined by some background metric g, and E"^^ is a non-degenerate symmetric 
matrix (a fiber metric). As usual, we assume that a summation over repeated 



M = Ix'Z, 



(2.1) 




(2.2) 



M 
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indices is performed. This metric (and its inverse E~^^^) can be used to naturally 
identify the bundle with its dual "V* (that is to raise and lower the field indices). 
The sections of the dual bundle are called currents and are represented locally by 
a set of functions, e.g. 

Ja=Eabv''. (2.3) 

We will also use the condensed De Witt notation, where the discrete index A and 
the spacetime point x are combined in one lower case Latin index i = {A,x). Then 
the components of a field ip are (i^') = {ifA{x)). There is a natural pairing between 
the bundles "V and "V* defined by 

{J,ip) = J,ip' = J dvoKx)JA{x)<p^{x) . (2.4) 

M 

It is assumed that a summation over repeated lower case Latin indices, i.e. a 
combined summation-integration, is performed. 

The set of all sections of the vector bundle is called the configuration space, 
which one assumes to be an infinite-dimensional manifold M. The fields (p' are the 
coordinates on this manifold, the variational derivative 6/6(p is a tangent vector, 
a small disturbance 6ip is a one-form and so on. If S{(p) is a scalar field on the 
configuration space, then its variational derivative 6S/6(p is a one-form on M 
defined by 



^S{<p + sh) 
as 



J«,L^^,.. (2.5, 

£=0 \6(p I 6(p' 



By using the functional differentiation one can define formally the concept of 
tangent space, the tangent vectors. Lie derivative, one-forms, metric, connection, 
geodesies and so on (for more details, see ^2M ). 



2.1 Non-Gauge Field Theories 

The dynamics of quantum field theory is determined by an action functional S {(f), 
which is a differentiable real-valued scalar field on the configuration space. The 
dynamical field configurations are defined as the field configurations satisfying 
the stationary action principle, i.e. they must satisfy the dynamical equations of 
motion 

^=0 (2.6) 

0(p 

with given boundary (and initial) conditions. The set of all dynamical field con- 
figurations, i.e. those that satisfy the dynamical equations of motion. Mo, is a 
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subspace of the configuration space called the dynamical subspace (or the mass- 
shell in the high-energy physics jargon). 

Quantum field theory is basically a theory of small disturbances on the dynam- 
ical subspace. Most of the problems of standard quantum field theory deal with 
scattering processes, which are described by the transition amplitudes between 
some well defined initial and final states in the remote past and the remote future. 
The collection of all these amplitudes is called the scattering matrix, or shortly 
S -matrix. 

Let us single out in the space-time two causally connected in- and out- re- 
gions, that lie in the past and in the future respectively relative to the region Q., 
which is of interest from the dynamical standpoint. Let |in) and |out) be some ini- 
tial and final states of the quantum field system in these regions. Let us consider 
the transition amplitude <out|in) and ask the question: how does this amplitude 
change under a variation of the interaction with a compact support in the region 
Q.. The answer to this question gives the Schwinger variational principle which 
states that 



where 6S is the corresponding change of the action. This principle gives a very 
powerful tool to study the transition amplitudes. The Schwinger variational prin- 
ciple can be called the quantization postulate, because all the information about 
quantum fields can be derived from it. 

Let us change the external conditions by adding a linear interaction with some 
external classical sources J in the dynamical region Q., i.e. 



The amplitude <out|in) becomes a functional of the sources that we denote by Z{J). 
The primary objects of interest in quantum field theory are the chronological mean 
values 



where T denotes the operator of chronological ordering that orders the (non- 
commuting) operators in order of their time variables from right to left. Of course, 
in the presence of the sources they become functionals of J. By using the Schwinger 
variational principle one can obtain the chronological mean values in terms of the 



5<out|in) = -<out|55' |in). 



fi 



(2.7) 



6S={J,<p). 



(2.8) 




(2.9) 
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functional derivatives of the functional Z{J), that is, 

z(y+77) = z(y)|i+(^)<77,^i(y)) + i(^)'<77,^2W 

+ Z;^(^)^«''(-^>'?'V--'7j. (2.10) 

n=3 ' 

In other words, the functional Z{J) is the generating functional for the chronolog- 
ical amplitudes *P„. 

Let us now define another functional W{ J) by 

Z = exp|^w). (2.11) 

Its functional derivatives define so called full connected Green functions, Qn"'^", 
(or the correlation functions) by 

1 °° 1 

W(7 + 7/) = W^(7) + (77,^1(7)) + - {n,&l{J)ri) + V -:^;] -'"(7)r/,, • • • 77/„ ■ (2. 12) 

z /z. 

n=3 

The functional = is called the background (or the mean) field, the operator 
^ = ^2 is called the full propagator. Then, it is easy to see that all chronological 
mean amplitudes can be expresssed in terms of connected Green functions. In 
particular, we have 

^1 = (p, (2.13) 
= <pjcl>^ + -^K (2.14) 

Thus, whilst Z(7) is the generating functional for chronological amplitudes, the 
functional W(J) is the generating functional for the connected Green functions. 
The Green functions satisfy the boundary conditions which are determined by the 
states |in) and |out). 

The mean field itself is a functional of the sources, (p = cf)(J). It is easy to see 
that the functional derivative of the mean field is equal to the full propagator, that 
is, 

^<P(J + sr]) =Qn. (2.15) 



de 



e=0 



In the non-gauge theories the full propagator Q, which plays the role of the (infinite- 
dimensional) Jacobian, is non-degenerate. Therefore, one can change variables 
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and consider ^ as independent variable and / = (as well as all other function- 
als) as a functional of ^. 

There are many different ways to show that there is a functional r(^) such that 

W = 2M. (2.16) 

This functional is defined by 

<out|in) = exp [r + {J, 0)]} , (2. 17) 

or by the functional Legendre transform 

ncl>) = W{J{(t>))-{J{(t>),(P). (2.18) 

This is the most important object in quantum field theory. It contains all the 
information about quantized fields. The functional expansion of this functional 
reads 

1 '^1 

r(0 +h) = T{cp) - {m,h) - - {h,g{m)h) + y -rr /i...,„(<^)/i'» • • • h^^ . (2. 19) 

n=3 

Therefore, the first variation of F gives the effective equations for the background 
fields 

^ = -J. (2.20) 
d(p 

These equations replace the classical equations of motion and describe the effec- 
tive dynamics of the background field with regard to all quantum corrections. That 
is why F is called the effective action. 

Furthermore, the second derivative of F(0) determines the full propagator 



G-l 

The higher derivatives, T^i-i,^, determine the so-called full vertex functions (also 
called strongly connected, or one-particle irreducible, functions). In other words, 
F(0) is the generating functional for the full vertex functions. The full vertex 
functions together with the full propagator determine the full connected Green 
functions and, therefore, all chronological amplitudes and, hence, the S -matrix. 



Ivan G. Avramidi: Effective Action in Quantum Gravity 



8 



Thus, the entire quantum field theory is summed up in the functional structure of 
the effective action. 

One can obtain a very useful formal representation for the effective action in 
terms of functional integrals (called also path integrals, or Feynman integrals). 
A functional integral is an integral over the (infinite-dimensional) configuration 
space At. Although a rigorous mathematical definition of functional integrals is 
absent, they can be used in perturbation theory of quantum field theory as an ef- 
fective tool, especially in gauge theories, for manipulating the whole series of 
perturbation theory. The point is that in perturbation theory one encounters only 
functional integrals of Gaussian type, which can be well defined effectively in 
terms of classical propagators and vertex functions. The Gaussian integrals do not 
depend much on the dimension and, therefore, (after a proper normalization) all 
formulas from the finite-dimensional case, like Fourier transform, integration by 
parts, delta-function, change of variables etc., are valid in the infinite-dimensional 
case as well. One has to note that functional integrals are formally divergent — if 
one tries to evaluate the integrals, one encounters meaningless divergent expres- 
sions. This difficulty can be overcome in the framework of the renormalization 
theory (in so-called renormalizable field theories). In non-renormalizable theo- 
ries (like quantum general relativity) this issue becomes the main difficulty of the 
theory. 

Integrating the Schwinger variational principle one can obtain the following 
functional integral: 



Here Dip represents the functional measure; however, it should not be taken too 
seriously — it will just provide a formal device for manipulations of Gaussian inte- 
grals. Correspondigly, for the effective action one obtains the functional equation 



M 

The only way to get numbers from this formal expression is to take advantage 
of the semi-classical approximation within a formal (asymptotic) expansion in 
powers of Planck constant h: 




(2.21) 



M 





00 




(2.23) 



k=\ 
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Next, we substitute this expansion in the functional equation for the effective ac- 
tion, shift the integration variable in the functional integral 

ip = (p+yfhh, (2.24) 

and expand the action S {(p) in functional Taylor series in quantum fields h 

00 

+ V-r/25,,-,..,„(0)/i'«---/i\ (2.25) 

where L is a (usually, partial differential) operator defined by the second variation 
of the action 

L = -^. (2.26) 

Notice that the operator L maps sections of the vector bundle 'V to sections of the 
dual bundle 'V*, that is, 

L : C~(^) ^ C'^(^*) . (2.27) 

In order to have a well defined operator which is self-adjoint with respect to the 
inner product on the bundle 'V we define another operator 

L:C°°('V)^C°°('V), (2.28) 

such that 

((p,Lh) = J dwol ip^EAsL^ch^ = J dwol (f^LAch^ = {(p,Lh) , (2.29) 

M M 

that is, 

EABL^c = LAch^. (2.30) 

Now, by expanding both sides of the functional equation for the effective ac- 
tion in powers on h and equating the coefficients of equal powers of h, we get 
the recurrence relations that uniquely define all coefficients Tg^y The measure for- 
mally transforms as D(p = Dh . All functional integrals appearing in this expansion 
are Gaussian and can be calculated in terms of the functional determinant, DetL, 
of the operator L and the bare propagator G = L~^, i.e. the Green function of the 
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operator L with Feynman boundary conditions. More precisely, with the proper 
normalization of the measure one can define 

exp(-^(/i,L/i)j = (DetL)"i/2^ (2.31) 

M 

JmQxp(^-^ih,Lhfjh''---h'^"'^' = 0, (2.32) 

M 

f DhQxp(-Uh,Lh)]h''---h'^'" = P'^]^' (DetL)~^/^G^'i'^ • • • G'^'"-i'^»'^ , 



M 



(2.33) 



where parenthesis denote the complete symmetrization over all indices included. 
Of course, the Green functions of the operators L and L are related by 

G^B(x,y) = G''^(x,y)EcB(y) ■ (2.34) 

In particular, the one-loop effective action is determined by the functional de- 
terminant of the operator L 

r(i) = -^logDetL, (2.35) 
and the two-loop effective action is given by 

r(2) = --S ^ijk0^G^^ - —S jjkG''^G^'"G'^"S jrnn ■ (2.36) 

Strictly speaking, the Gaussian integrals are well defined for elliptic partial 
differential operators in terms of the functional determinants and their Green func- 
tions. Although the Gaussian integrals of quantum field theory are determined by 
hyperbolic partial differential operators with Feynman boundary conditions they 
can be well defined by means of the analytic continuation from the Euclidean sec- 
tor of the theory where the operators become elliptic. This is done by so-called 
Wick rotation — one replaces the real time coordinate by a purely imaginary one 
jc^ ir and singles out the imaginary factor also from the action S iS and the 
effective action F /F. Then the metric of the spacetime manifold becomes posi- 
tive definite and the classical action in all 'nice' field theories becomes a positive- 
definite functional. Then the fast oscillating Gaussian functional integrals become 
exponentially decreasing and can be given a rigorous mathematical meaning. 
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2.2 Gauge Field Theories 

Let us try to apply the formalism described above to a gauge field theory. A char- 
acteristic feature of a gauge field theory is the fact that the dynamical equations 

6S 

— = (2.37) 

0(fi 

are not independent — there are certain identities, called Nother identities, be- 
tween them. This means that there are some nowhere vanishing vector fields 

%, = R\-^, (2.38) 

on the configuration space M that annihilate the action, 

RaS=Q, (2.39) 

and, hence, define invariance flows on M. The transformations of the fields 

6^<p'=R\e' (2.40) 

are called the invariance transformations and Rq, are called the generators of in- 
variance transformations. The infinitesimal parameters of these transformations ^ 
are sections of another vector bundle (usually the tangent bundle TG of a com- 
pact Lie group G) that are respresented locally by a set of functions {^") - (^"(^)) , 
a = 1, . . . ,dimG, over spacetime with compact support. To distinguish between the 
components of the gauge fields and the components of the gauge parameters we 
introduce lower case Latin indices from the beginning of the alphabet; the Greek 
indices from the beginning of the alphabet are used as condensed labels a = (a,x) 
that include the spacetime point. 

We assume that the vector fields Rq. are linearly independent and complete, 
which means that they form a complete basis in the tangent space of the invariant 
subspace of configuration space. The vector fields Ra form the gauge algebra. We 
restrict ourselves to the simplest case when the gauge algebra is the Lie algebra of 
an infinite-dimensional gauge Lie group ^. This is the case in Yang-Mills theory 
and gravity. Then the flow vectors Rq, decompose the configuration space into the 
invariants subspaces of M (called the orbits) consisting of the points connected by 
the gauge transformations. The space of orbits is then MIQ. The linear indepen- 
dence of the vectors Rq, at each point implies that each orbit is a copy of the group 
manifold. One can show that the vector fields Rq, are tangent to the dynamical 
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subspace Mq, which means that the orbits do not intersect Mq and the invariance 
flow maps the dynamical subspace Mq into itself. Since all field configurations 
connected by a gauge transformation, i.e., the points on an orbit, are physically 
equivalent, the physical dynamical variables are the classes of gauge equivalent 
field configurations, i.e., the orbits. The physical configuration space is, hence, 
the space of orbits MIQ. In other words the physical observables must be the 
invariants of the gauge group. 

To quantize a gauge theory by means of the functional integral, we consider the 
in- and out- regions, define some |in) and |out) states in these regions and study 
the amplitude (out|in). Since all field configurations along an orbit are physically 
equivalent we have to integrate over the orbit space MIQ. To deal with such situ- 
ations one has to choose a representative field in each orbit. This can be done by 
choosing special coordinates {I^{ip),x"{''P)) on the configuration space M, where 
label the orbits and the points in the orbit. Computing the Jacobian of the 
field transformation and introducing a delta functional 6{x - we can fix the co- 
ordinates on the orbits and obtain the measure on the orbit space MIQ 

m = DcpDctFi<p)6(x(cp) - 0, (2.41) 

where 

F''a=^a/ (2.42) 

is a non-degenerate operator. Thus we obtain a functional integral for the transi- 
tion amplitude 

<out|in) = j£>^DetF(^)(J(;kr(^)-^)exp{^5(^. (2.43) 
M 

Now one can go further and integrate this equation over parameters ^ with a 
Gaussian measure determined by a symmetric nondegenerate matrix 7 = {jap), 
which most naturally can be chosen as the metric on the orbit (gauge group met- 
ric). As a result we get 



<out|in) = ^ D(p{T>etyfl^ DetF(^)exp|^ 
M 



(2.44) 



The functional equation for the efi'ective action takes the form 

exp|^r((^)j = Jl)^ (Det7(0))i/2 DetF(^) (2.45) 



M 



Ivan G. Avramidi: Effective Action in Quantum Gravity 



13 



xexp^l 



The determinants of the operators F and y are usually represented as a result of 

the integration over some auxiliary Grassmanian variables, so called ghost fields. 

This equation can be used to construct the semi-classical perturbation theory 
in powers of the Planck constant (loop expansion), which gives the effective action 
in terms of the bare propagators and the vertex functions. In particular, one finds 
the one-loop eff'ective action 

1 .1 1 

Td) = -— logDet L+-logDetF + — logDety , (2.46) 
2i I 2i 



where L is an operator defined by 

= -(h,Lh). (2.47) 

I e=0 

In De Witt notation it reads 



j2 r 1 

-^<S((p + Eh) + - {x((p + eh), yx{(p + eh)) 



L^^. = (2.48) 

where 

6^S (Jy" 5x^ 

Lij = : : - -^JaB^ . (2.49) 

2.3 Quantum General Relativity 

Einstein's theory of general relativity is an example of a gauge theory with the 
gauge group ^ being the group of all diffeomorphisms of the spacetime manifold 
M and the configuration space M being the space of all pseudo-Riemannian met- 
rics on M. The physical configuration space M/Q of all orbits of the gauge group 
is then the space of all geometries on the spacetime. 

The gravitational field can be parametrized by the metric tensor of the space- 
time 

(fi' = g,_ivix) , i = {jiv,x). (2.50) 
An invariant fiber metric is defined by 

E^^vafi ^ gn{a^)v _ ^^y^ap ^ (2.5 1) 

where xi^llnisdL real parameter. The inverse metric is then 

-1 ^ 

^Hvap = 8fi(ag/3)v - -;—[Siivgap ■ (2.52) 
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The parameters of gauge transformations are the components of the vector of 
the infinitesimal diffeomorphism, 

e=e{x), ^i^(M,x). (2.53) 

An invariant metric in the gauge group can be chosen to be just a background 
metric g^y. 

The local generators of the gauge transformations in this parametrization are 
defined by their action as follows 

= 2V(^^,), i = (MV,x), (2.54) 

JiR'a = -2V^J^a , a = (a,x). (2.55) 
The Hilbert-Einstein action of general relativity has the form 

dxg^'^(R-2A) , (2.56) 

M 

where R is the scalar curvature, = \6nG is the Einstein coupling constant, G 
is the Newtonian gravitational constant and A is the cosmological constant. Here 
we neglect the boundary term for simplicity; it will not affect our calculations. 
The first variation of the action gives the classical equations of motion 



?-'»^ = -;^(i?'"'-ij;'"'fi + Ar). (2.57) 

which satisfy, of course, the Nother identities 

V^|i?^^-ig^^i? + Ag^^j = 0. (2.58) 

Here R^y is the Ricci tensor defined in terms of the Riemann tensor by R^y = R"j_iav 
The second variation of the action defines a second-order partial differential 
operator by 

= (2.59) 

OgnvOgap 

where 
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+ (g^ - j - 2A)| (2.60) 

Here, of course, A = g^^V^Vy denotes the Laplacian. 
Next, we choose the De Witt gauge condition 

= -^"^'^^v^v = -(g'^'v''^ -^g'^'v'^) V • (2.61) 

The ghost operator in this gauge is a second-order differential operator defined by 

F^y = -lE^^l^yVa'^p = -(5^^ A + (2x - l)V^Vv - < . (2.62) 

For this operator to be non-singular, the gauge parameter should satisfy the con- 
dition x 1. 

For the graviton operator L to be non-degenerate it is necessary to choose the 
operator y as a zero order differential operator defined by 

7nv = ^gfiv , (2.63) 

where or is a real parameter. Thus we obtain a two-parameter class of gauges 
involving two arbitrary parameters, x and a. 
The graviton operator L now reads 

^^^'"^ = ^|-(g'^V^-(l+2axVV> 

-(1 + 2a;<)g^^V^«V^) - (1 + 2ax)g«^V^V^) + 2(1 + a)V<^g^)(«V^^ 

+ (g^ - V^) (R - 2A)} (2.64) 
The most convenient choice is the so-called minimal gauge 

K = ^, a = -l. (2.65) 
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In this gauge the non-diagonal derivatives in both the graviton operator and the 
ghost operator vanish 



F^'y = -6^yA-R^'y. (2.67) 

Finally, we define the graviton operator in the canonical Laplace-type form, L, 
by factoring out the configuration space metric (in the minimal gauge x = l/2) 

L^/P = lk%lp,LP--^. (2.68) 



We obtain 
where 



One can show that the contribution of the determinant of the operator y can be 

neglected (more precisely, it can be absorbed in the measure of the path integral) 
since it is of zero order. Thus, with this choice of gauge parameters the one-loop 
effective action of quantum general relativity is given by 

1 . 1 

Tn) = -— logDetL+-logDetF. (2.71) 

2i I 

Therefore, in order to compute the eff"ective action we need to compute the de- 
terminants of Laplace type partial differential operators acting on symmetric two- 
tensors and vectors. 



3 Heat Kernel Method 



As we described in the previous section the effective action in quantum field theory 
can be computed within the semi-classical perturbation theory. It is determined 
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by the functional determinants of second-order hyperbolic partial differential op- 
erators with Feynman boundary conditions and the higher-loop approximations 
are determined in terms of the Feynman propagators and the classical vertex func- 
tions. As we noted above these expressions are purely formal and need to be 
regularized and renormalized, which can be done in a consistent way in renor- 
malizable field theories. One should stress, of course, that many physically in- 
teresting theories (including Einstein's general relativity) are perturbatively non- 
renormalizable. Since we only need Feynman propagators we can do the Wick 
rotation and consider instead of hyperbolic operators the elliptic ones. The Green 
functions of elliptic operators and their functional determinants can be expressed 
in terms of the heat kernel. That is why we concentrate below on the calculation 
of the heat kernel. 

The heat kernel is one of the most powerful tools in mathematical physics and 
geometric analysis (see, for example, the books [STl [TSl [TOl l28l and reviews |fT9l 
|9l[ITl[32l[14l|). The short-time asymptotic expansion of the trace of the heat kernel 
determines the spectral asymptotics of the differential operator. The coefficients 
of this asymptotic expansion, called the heat invariants, are extensively used in 
geometric analysis, in particular, in spectral geometry and index theorems proofs 

mm. 

The gauge invariance (or covariance) in quantum gauge field theory and quan- 
tum gravity is of fundamental importance. That is why, manifestly covariant meth- 
ods present inestimable advantage. A manifestly covariant calculus is such that 
every step is expressed in terms of geometric objects; it does not have some in- 
termediate non-covariant steps that lead to an invariant result. Below we describe 
a manifestly covariant method for calculation of the heat kernel following mainly 
our papers [l2l[3l|9l[IOl[IIl. 

3.1 Laplace Type Operators 

Let iM,g) be a smooth compact Riemannian manifold of dimension n without 
boundary, equipped with a positive definite Riemannian metric g. We assume 
that it is complete simply connected orientable and spin. We denote the local 
coordinates on M by x^, with Greek indices running over \,...,n. Let be a 
local orthonormal frame defining a basis for the tangent space T^M. We denote 
the frame indices by low case Latin indices from the beginning of the alphabet, 
which also run over l,...,n. The frame indices are raised and lowered by the 
metric 6ah- Let e"^ be the matrix inverse to ea^, defining the dual basis in the 
cotangent space r*M. As usual, the orthonormal frame, e"^ and e^^, will be 
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used to transform the coordinate (Greek) indices to the orthonormal (Latin) in- 
dices. The Riemannian volume element is defined as usual by Jvol = dxg^^^ , 
where g = det^^y = (detea^)^ . The spin connection tj"''^ is defined in terms of the 
covariant derivatives of the orthonormal frame with the Levi-Civita connection. 
The curvature of the spin connection defines the Riemann tensor, R^bnv^ the Ricci 
tensor, R^y = R^^av, and the scalar curvature, R = R^^, as usual. 

Let 7~ be a spin-tensor bundle realizing a representation S of the spin group 
Spin(n), the double covering of the group SO{n), with the fiber A. Let "Lab be 
the generators of the orthogonal algebra SO{n), the Lie algebra of the orthogonal 
group S 0(n). The spin connection induces a connection on the bundle T defining 
the covariant derivative of smooth sections (p of the bundle T by 

Vr^ = (5^ + ^6>«^S«,j^. (3.1) 

The commutator of covariant derivatives defines the curvature of this connection 
via 

[vr,vr]^=ii?«V2„,^. (3.2) 

The covariant derivative along the frame vectors is defined by = eJ^V ^. For 
example, with our notation, V oStTcd = ^ct^b ^"^(P'^ yTafi- The metric 5ab 
induces a positive definite fiber metric on tensor bundles. 

Let Gym be a compact Lie group (called a gauge group). It naturally defines 
the principal fiber bundle over the manifold M with the structure group Gym- We 
consider a representation of the structure group Gym and the associated vector 
bundle through this representation with the same structure group Gym whose typ- 
ical fiber is a ^-dimensional vector space W. Then for any spin-tensor bundle T 
we define the twisted spin- tensor bundle 'V via the twisted product of the bundles 

and T. The fiber of the bundle "V is V = A (g) W so that the sections of the 
bundle "V are represented locally by ^-tuples of spin-tensors. 

Let be a connection one form on the bundle "W (called Yang-Mills or 
gauge connection) taking values in the Lie algebra Qym of the gauge group Gym- 
Then the total connection on the bundle is defined by 

Vf,(p = (d^,+:R^)ip, (3.3) 

where 

= ^Co'^^^llab^Iw + lA^^l^, (3.4) 
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and the total curvature H of the bundle 'V is defined by 



(3.5) 



where 



'■/j.v 



(3.6) 



and 




(3.7) 



is the curvature of the Yang-Mills connection. 

We also consider the bundle End('y) of endomorphisms of the bundle 'V. The 
covariant derivative of sections of this bundle is defined by 



We assume that the vector bundle "V is equipped with a Hermitian metric. This 
naturally identifies the dual vector bundle "V* with "V. We assume that the con- 
nection V is compatible with the Hermitian metric on the vector bundle 'V. The 
connection is given its unique natural extension to bundles in the tensor algebra 
over and "V*. In fact, using the Levi-Civita connection of the metric g together 
with the connection on the bundle "V, we naturally obtain connections on all bun- 
dles in the tensor algebra over "Y, 'V*, TM and T*M; the resulting connection will 
usually be denoted just by V. It is usually clear which bundle's connection is being 
referred to, from the nature of the section being acted upon. 

We denote by C^CV) the space of smooth sections of the bundle 'V. The fiber 
inner product on the bundle 'V defines a natural inner product and the -trace 
Tr using the invariant Riemannian measure on the manifold M. The completion 
of C°°('V) in this norm defines the Hilbert space L^CV) of square integrable sec- 
tions. Let V* be the formal adjoint to V defined using the Riemannian metric 
and the Hermitian structure on 'V and let Q be a smooth Hermitian section of the 
endomorphism bundle End('y). 

A Laplace type operator L : C°°(V) C°°(V) is a partial differential operator 
of the form 



(3.8) 



and the commutator of covariant derivatives is equal to 



[V.,v,](2 = [^.v,!2]. 



(3.9) 



L = TV + Q = -A + Q. 



(3.10) 
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In local coordinates the Laplacian is defined by 

A =g^%V, =^-1/2(5^ +j?l^)^l/2^/'V(5^+J?l^). (3.11) 

and, therefore, 

= -g^'d^dy-la^d^ + q, (3.12) 

where 

a^'=g^'My + ^g-'/^dyig'/^g'^') (3.13) 

q = Q-g^''Ji^Jly-g-'^%{g'/V~^y). (3.14) 
Thus, a Laplace type operator is constructed from the following three pieces 
of geometric data: i) a Riemannian metric g on M, which determines the second- 
order part, ii) a connection 1-form Jl on the vector bundle "V, which determines 
the first-order part, iii) an endomorphism Q of the vector bundle "V, which deter- 
mines the zeroth order part. It is worth noting that every second-order differential 
operator with a scalar leading symbol given by the metric tensor is of Laplace 
type and can be put in this form by choosing the appropriate connection and the 
endomorphism Q. 

It is easy to show that the Laplacian, A, and, therefore, the operator L, is an 
elliptic symmetric partial differential operator satisfying 

(L^,(A) = (<^,L(A), (3.15) 

with a positive principal symbol. Moreover, the operator L is essentially self- 
adjoint, i.e., it has a unique self-adjoint extension. We will not be very careful 
about distinguishing between the operator L and its closure, and will simply say 
that the operator L is elliptic and self-adjoint. 
It is well known [27J that: 

i) the operator L has a discrete real spectrum, {A„}'^^^, bounded from below: 

Ao < Ai < A2 < ■■■ < An < ■■■ (3.16) 
with some real constant Aq, 

ii) the eigenvalues grow as ^ 00 as /lyt ~ Ck^^", where n = dimM, 

iii) all eigenspaces of the operator L are finite-dimensional, and 

iv) the eigenvectors, of the operator L, are smooth sections of the vec- 
tor bundle that form a complete orthonormal basis in L^CV). 
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3.2 Spectral Functions 

The spectrum of the operator L can be described by cerain spectral invariants, 
called spectral functions. First of all, we define the heat trace: 

00 

where each eigenvalue is counted with multiplicities. The heat trace is well de- 
fined for real positive t. Notice that it can be analytically continued to an analytic 
function of t in the right half -plane (for Re ? > 0). 

The heat trace determines other spectral functions by integral transforms: the 
distribution function (also called counting function), defined as the number of 
eigenvalues below the level A, 

E+ioo 

\ r dt 

N{X) = Yj(^{A-Xn) = — \ je'^@(t), (3.18) 
where e is a positive constant, the density function, 

p(^) = 25(A-^) = — j dte'^mX (3.19) 
and the zeta-f unction, 

00 

°° 1 1 r 

^(^-^)^ T{s)j 

where /I is a large negative constant such that Re/l < /lo and 5 is a complex param- 
eter with Re 5 > nil. 

In principle, if known exactly, they determine the spectrum. Of course, this is 
not valid for asymptotic expansions of the spectral functions. There are examples 
of operators that have the same asymptotic series of the spectral functions but 
difi'erent spectrum. 

The zeta function enables one to define, in particular, the zeta-regularized de- 
terminant of the operator (L - A), 

C(0,A)^-^as,A.) =-logDet(L-i), (3.21) 

OS s^O 

which determines the one-loop eff'ective action in quantum field theory. 
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3.3 Heat Kernel 

For ? > the operators 

U{t) = exp(-?L) (3.22) 

form a semi-group of bounded operators on L^CV), so called heat semi-group. 
The kernel of this operator is defined by 

CO 

U(t\x,x') = J]e-''^''<p„(x)^<pl(x'), (3.23) 

n=l 

where each eigenvalue is counted with multiplicities. It is a section of the external 

tensor product of vector bundles 'V El "V* over MxM, which can also be regarded 
as an endomorphism from the fiber of 'V over x' to the fiber of 'V over x. This 
kernel satisfies the heat equation 

(dt + L)U(t) = (3.24) 

with the initial condition 

U(0^\x,x') = 6(x,x') (3.25) 

and is called the heat kernel. 

Moreover, the heat semigroup U(t) is a trace-class operator with a well defined 
-trace, 

Trexp(-?L) = J dwoltTvU'^''Ht). (3.26) 

M 

Hereafter try denotes the fiber trace and the label 'diag' means the diagonal value 
of a two-point quantity, e.g. 

U'^'^Ht\x) = U(t\x,x') . (3.27) 

x=x' 

It is easy to see that the heat trace defined above is equal to the trace of the heat 
semigroup, that is, 

0(0 = Tr exp(-?L) . (3.28) 



3.4 Asymptotic Expansion of the Heat Kernel 



In the following we are going to study the heat kernel only locally, i.e. in the 
neighbourhood of the diagonal of M x M, when the points x and x' are close to 
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each other. The exposition will follow mainly our papers |l3l[T0l|9l[IIl. We will 
keep a point of the manifold fixed and consider a small geodesic ball, i.e. a 
small neighbourhood of the point x': B^ix') = [x e M\r{x,}^) < s), r{x,x') being 
the geodesic distance between the points x and x' . We will take the radius of the 
ball sufficiently small, so that each point x of the ball of this neighbourhood can 
be connected by a unique geodesic with the point x' . This can be always done if 
the size of the ball is smaller than the injectivity radius of the manifold, s < r^y 

Let cr{x,x') be the geodetic interval, also called world function, defined as one 
half the square of the length of the geodesic connecting the points x and x' 

(rix,x') = ^r\x,x'). (3.29) 

The first derivatives of this function with respect to x and x' define tangent vec- 
tor fields to the geodesic at the points x and x' respectively pointing in opposite 
directions 

= g^^VyCT, (3.30) 
u^' = g^'^'v'^,cr, (3.31) 

and the determinant of the mixed second derivatives defines a so called Van Vleck- 
Morette determinant 

A{x,x')=g-^/\x)det[-Vi,Y^,(T{x,/)]g-^'\x'). (3.32) 

This object should not be confused with the Laplacian, which is also denoted by 
A. 

Let, finally, 'P{x,x') denote the parallel transport operator of sections of the 
vector bundle 'V along the geodesic from the point x' to the point x. It is a section 
of the external tensor product of the vector bundle "V El "V* over M x M, or, in 
other words, it is an endomorphism from the fiber of "V over x' to the fiber of 
over X. Here and everywhere below the coordinate indices of the tangent space at 
the point x' are denoted by primed Greek letters. They are raised and lowered by 
the metric tensor g^'v'(x') at the point x' . The derivatives with respect to x' will 
be denoted by primed Greek indices as well. 

We extend the local orthonormal frame ea^ (x') at the point x' to a local or- 
thonormal frame ej^ix) at the point x by parallel transport. The parameters of the 
geodesic connecting the points x and x' , namely the unit tangent vector at the point 
x' and the length of the geodesic, (or, equivalently, the tangent vector at the point 
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x' with the norm equal to the length of the geodesic), provide normal coordinate 
system for 5e(jc'). Now, let us define the following geometric parameters 

y- = e%u'' = -e%>u^'' , (3.33) 

so that 

M^=e// and M^' = -e/y. (3.34) 

Notice that = dX x = x' . The geometric parameters are nothing but the 
normal coordinates. 

Near the diagonal of M x M all these two-point functions are smooth single- 
valued functions of the coordinates of the points x and x' . Let us note from the 
beginning that we will construct the heat kernel in form of covariant Taylor series 
in coordinates. In the smooth case these series do not necessarily converge. How- 
ever, if one assumes additionally that the two-point funtions are analytic, then the 
Taylor series converge in a sufficiently small neighborhood of the diagonal. 

Further, one can easily prove that the function 

Uo{t\x,x') = (47rO~"^^Ai/2(x,x')exp|-^o-(x,x')j^(x,x') (3.35) 
satisfies the initial condition 

Uo(0^\x,x') = 6(x,x'). (3.36) 

Moreover, locally it also satisfies the heat equation in the free case, when the Rie- 
maimian curvature of the manifold, Riem, the curvature of the bundle coimection, 
and the endomorphism Q vanish: 

Riem = :R = 2 = 0. (3.37) 

Therefore, Uo(t\x, x') is the exact heat kernel for a pure Laplacian in flat Euclidean 
space with a flat trivial bundle connection and without the endomorphism Q. 

3.4.1 Transport Function 

This function gives a good framework for the approximate solution in the general 
case. Namely, by factorizing out this free factor we get an ansatz 

Ui,.X> - (4«)-A-(.,/,exp(-i.(„',)n.,/)n*,.',. (3.38, 
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The function 0.{t\x,x'), called the transport function, is a section of the endo- 
morphism vector bundle End(y) over the point x' . Using the definition of the 
functions cr{x,x'), A(;c,x') and 'P{x,x') it is not difficult to find that the transport 
function satisfies a transport equation 



Q(0 = 0, (3.39) 



where D is the radial vector field, i.e. operator of differentiation along the geodesic, 
defined by 

D = M^V^, (3.40) 
and L is a second-order differential operator defined by 

L = 9-^^-^>^L^^'^r. (3.41) 

The initial condition for the transport function is obviously 

Q.{t\xJ) = \, (3.42) 

where is the identity endomorphism of the vector bundle over x' . 

It is obvious that if we replace the operator L by (L-/1), with Rc/l < /Iq, then 
the heat kernel and the transport function are simply multiplied by e''^, i.e. the 
transport function for the operator (L - /I) is e^'^Oit). Further, for X< the oper- 
ator (L-/1) becomes a positive operator. Therefore, the function e^'^Vlit) satisfies 
the following asymptotic conditions 

lim t"d^\e'^a{tS\ =0 for i < ii, o- > 0, > 0. (3.43) 

In other words, as ? ^ oo the function e^'^£i{f) and all its derivatives decreases 
faster than any power of t, actually it decreases exponentialy, and as ? ^ the 
product of ^'"^^(0 with any positive power of t vanishes. 

Hereafter we fix /I < /Iq, so that (L - /I) is a positive operator. Now, let us 
consider a slightly modified version of the Mellin transform of the function e^'^^iif) 
introduced in [3J 

oo 

^^^-^^ = r ^''"'^^''^^(0- (3.44) 

r(-^) J 



Note that for fixed X this is a Mellin transform of e^'^Oit) and for a fixed q this 
is a Laplace transform of the function ?~^~^Q(?). The integral (13.441) converges 
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for Re^ < 0. By integrating by parts times and using the asymptotic conditions 
(13.431 ) we also get 

oo 

= Y{- \n) / ^-^'^^ • (3.45) 



This integral converges for Re^ < A'^ - 1. Using this representation one can prove 
that [3] the function bq{A) is an entire function of q (analytic everywhere) satisfy- 
ing the the asymptotic condition 

lim Yi-q + N)bq{A) = 0, for any > 0. (3.46) 

\q\—^oo, Req<N 

Moreover, the values of the function bq(A) at the integer positive points q = k axe 
given by 

bkW = {-dtf [e'^Q{t)] = E Q^" ' (3-^^^ 

where 

a^ = (-5/Q(0 , (3.48) 

By inverting the Mellin transform we obtain a new ansatz for the transport 
function and, hence, for the heat kernel 



2ni / 



C+lOO 



Q(0 = — dqe-'^tlT{-q)bM) (3.49) 



where c < and 'Re.A< Aq. Clearly, since the left-hand side of this equation does 
not depend on A, neither does the right hand side. Thus, A serves as an auxiliary 
parameter that regularizes the behavior at ? ^ oo. If we invert instead the Laplace 
transform, we obtain another representation 

7+;'oo 

Q(0 = — r dAe-'^t'i^^Y{-q)bq{A) (3.50) 
Ini J 

-y—ioo 

where y < Aq and Req <0. 

Substituting this ansatz into the transport equation we get a functional equation 
for the function bq 

l + ^Dyq(A) = iL-A)bq-i{A). (3.51) 
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The initial condition for the transport function is translated into 

bQ{A) = Iv. (3.52) 

Thus, we have reduced the problem of solving the heat equation to the fol- 
lowing problem: one has to find an entire function of q, bq{A\x,x'), that satisfies 
the functional equation (13.511) with the initial condition (13.521 ) and the asymptotic 
condition (13.461 ). 

Although the variables q and A seem to be independent they are very closely 
related to each other. In particular, by diff"erentiating with respect to A we obtain 
an important result 

— bq{A) = -qbq-,{A). (3.53) 

Also, by diff"erentiating the eq. (13.511) with respect to q one obtains another recur- 
sion 

|l + -^D^b'^iX) = Lb'^^M) + ^^h{A\ (3.54) 

where 

b'q{A) = -^^bqiAX (3.55) 

which enables one to compute the derivatives of the function bq{A) at positive in- 
teger points if one fixes its value b'^^A). This tims out to be useful when computing 
the determinant of the operator {L- A). 

Moreover, one can actually manifest the dependence of bq{A) on A. It is not 
difficult to prove that [3J the integral 



Y{ -p)Y{p- q)^ 

n-q) 

ci-ioo 



^^^-^^ = 2^ / (3.56) 



with Re^ < ci < 0, satisfies the equation (13.511) if ap satisfies this equation for 
i = 0, i.e. 

l + ^Dyq = Laq-i. (3.57) 
with the initial condition 

ao = lv. (3.58) 

For integer q = k = 1,2,... the functional equation (13.571) becomes a recursion 
system that, together with the initial condition (|3.58l) . determines all coeflicients 
ak- 
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Now, from eq. (13.561) we also obtain the asymptotic expansion of bq{A) as 
A -oo, 



For integer q this coincides with (13.471) . 

The function bq{A) turns out to be extremely useful in computing the heat 
kernel, the resolvent kernel, the zeta-function and the determinant of the operator 
L. It contains the same information about the operator L as the heat kernel. In 
some cases the function bq{A) can be constructed just by analytical continuation 
from the integer positive values bk 01 • 

3.4.2 Asymptotic Expansion of the Transport Function 

Now we are going to do the usual trick, namely, to move the contour of integration 
over q to the right. Due to the presence of the gamma function r(-^) the integrand 
has simple poles at the non-negative integer points ^ = 0, 1, 2 . . ., which contribute 
to the integral while moving the contour. So, we get 



with Cm is a constant satisfying the condition N - I < c^ < N. As t ^ the rest 
term R^it) behaves like 0(t^), so we obtain an asymptotic expansion as ? ^ 




(3.59) 




(3.60) 



where 




(3.61) 




(3.62) 



Using our ansatz (13.381 ) we find immediately the heat trace 




(3.63) 



c-ioo 
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where 



Bq{A) = ^ dvo\\.xvb^q''\X). (3.64) 



M 

The heat trace has an analogous asymptotic expansion as ? ^ 



0(0 - (4nt)-"'^e-'' J] ^BkW = J] , (3.65) 

k=0 ' k=0 



fc=U 

where 



Ayt = J" Jvol trva^^^^. (3.66) 
M 

This is the famous Minakshisundaram-Pleijel asymptotic expansion. The physi- 
cists call it the Schwinger-De Witt expansion [14J. Its coefficients are also 
called sometimes Hadamard-Minakshisundaram-De Witt-Seeley (HMDS) coeffi- 
cients. This expansion is of great importance in differential geometry, spectral 
geometry, quantum field theory and other areas of mathematical physics, such as 
theory of Huygens' principle, heat kernel proofs of the index theorems, Korteveg- 
De Vries hierarchy, Brownian motion etc. 

One should stress, however, that this series does not converge, in general. In 
that sense our ansatz (13.491 ) or (13.601) in form of a Mellin transform of an entire 
function is much better since it is exact and gives an explicit formula for the rest 
term. 



3.5 Zeta Function and Determinant 

Let us apply our ansatz for computation of the complex power of the operator 
(L - A) (with /I < /lo so that the operator (L - A) is positive) defined by 

CO 

GsiA) = (L-A)-' = -^A dtt'-' e'^ U(t). (3.67) 

r(5) J 



Using our ansatz for the heat kernel one can obtain ||3l 

c+ioo 

GAA)-i4.r-'^A'nr± f . /<-^>^'-^-""^^> (fp^,a) (3.68) 
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where c < -Re p + njl. 

Outside the diagonal, i.e. for cr 0, this integral converges for any s and 
defines an entire function of s. The integrand in this formula is a meromorphic 
function of q with some simple and maybe some double poles. If we move the 
contour of integration to the right, we get contributions from the simple poles 
in form of powers of cr and a logarithmic part due to the double poles (if any). 
This gives the complete structure of diagonal singularities of Gs{x,x'). Thus the 
function bq{X) turns out to be very useful to study the diagonal singularities. 

Now, let us consider the diagonal limit of G^. By taking the limit cr ^ we 
obtain a very simple formula in terms of the function bq 

Gf'HX) = (4;r)-/2lI^lZp^^diai_^(^). (3.69) 

This gives automatically the zeta-function 

Y(s — n/2) 

as,A) = (4;r)-"/" ' ' Bn/2-s(A). (3.70) 

r{s) 

Herefrom we see that both Gf^^(A) and ^(s,A) are meromorphic functions 
of s with simple poles at the points s = [n/2] + 1/2- k, (k = 0,1,2,...) and s = 
1,2, . . . , [n/2]. In particular, the zeta-function is analytic at the origin. Its value at 
the origin is given by 

for odd n , 

(4^) 77-^mT^«/2(^) for even n. 
r(n/2 + 1) 

This gives the regularized number of all modes of the operator L. 

Moreover, the derivative of the zeta-function at the origin is also well defined. 
As we already mentioned above it determines the regularized determinant of the 
operator (L - A) 

,^W_1^(n+l)/2 

logDet(L-^) = -(471)-"'^ ; ; B„/2(A) (3.72) 

r(n/2 + l) 

for odd n, and 

logDet(L -A) = (4;r)-"/2_±^!^l_ {b'^^^(A) - [¥(n/2 + 1) + C]5„/2(A)j (3.73) 
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for even n. Here = {dldz)\ogY{z) is the psi-function, C = -^(1) = 0.577 . . . 
is the Euler constant, and 

(3.74) 

q=nl2 

4 Green Function 

In this section we closely follow our paper fS^ . Let /I be a sufficiently large nega- 
tive parameter, such that A< Aq and, therefore, (L - /I) be a positive operator. The 
Green function of the operator (L - A) reads 

G{A\x,x') = y -<fn(.x)^^:{x'). (4.1) 

It is not difficult to see that the Green function can be represented as the Laplace 
transform of the heat kernel 

00 

G(A) = J dte'^U{t). (4.2) 



Using our ansatz for the heat kernel we obtain 

c+z'oo 

G(^) = (47r)-"/2^i/¥^ J dqY{-q)Y{-q-\+nl2){^Y " W (4-3) 

C— iOO 

where c < n/2 - 1. 

This ansatz is especially useful for studying the singularities of the Green func- 
tion, or more general, for constructing the Green function as a power series in a. 
The integrand in (|4.3I) is a meromorphic function with poles at the points q = k 
and q = k- \ + n/2, where (fc = 0, 1 , 2, . . .). Here one has to distinguish between 
odd and even dimensions. In odd dimensions, the poles are at the points q = k and 
^ = ^-1- [n/2] -1/2 and are simple, whereas in even dimension there are simple 
poles at ^ = 0, 1,2, . . . ,n/2 -2 and double poles at the points q = k + n/2- I. 

Moving the contour of integration in (14.31) to the right one can obtain an ex- 
pansion of the Green function in powers of cr (Hadamard series). Generally, we 
obtain 

G(A) = G'^'^HA) + G"°"-^"^i(^) + G'^H^) . (4.4) 



d_ 

dq 
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Here G^™^{X) is the singular part which is polynomial in the inverse powers of 

[{n+\)l2\-2 ^. jryyill-k-l 

G''^HX) = m-"'^^"^r g L_Lr(n/2-fc-l)|-j bk{A\ (4.5) 

Let us fix an integer A'^ such that N > (n- l)/2. 
For the rest we get in odd dimensions 



N-(n+l)/2 



Cf^ + lOO 

+(4n)-"'^A"^'P^. J dql^^J^ " r(-q)Yi-q-l+n/2)bqiA), 

Cff-ioo 

(4.6) 

where - 1 < c^ < N - 1/2. Thus, by putting N ^ oo we recover the Hadamard 
power series in cr for odd dimension n 

G-"\A) - (-1)("-1)/2(4;,)-/2a1/2^ y z.,^^(i) 

(4.7) 

G-S(^) ~ (-l)("-^l)/2(4;r)-"/2Ai/¥^ fc!r(fcTn/2) (f ) ^^-l-/2(^). (4.8) 

In even dimensions, the point is more subtle due to the presence of double 
poles. Moving the contour in (14.31) to the right and calculating the contribution of 
the residues at the simple and double poles we obtain 

G"°"-^""i(i) + G'"'=S(^) 



k=0 
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.(-,r'-(4.r'^A.'vt:^^(^)' 



^ k\Y{k + nl2)\2, 



X 



K-i+ni2^^) - [log)"' +^>{k+l) + ^>{k + nil)] bk_,^n{X) 



Cf^ + lOO 

+(4;r)-"/2Ai/2^2. J dg I^^J^'^"'\(-q)n-q- I +n/2)bq{A), 



(4.9) 



where is an arbitrary mass parameter introduced to preserve dimensions, A^- 1 < 
cn < N and ¥(z) = (d/dz)logr(z). If we let A^^ ^ oo we obtain the Hadamard 
expansion of the Green function for even dimension n>2 

^)l. kink^n/2) {^) ''-'^"''^'^ 

(4.10) 

°° 1 k 

G-H,> ~ (-,y.'-,4.r"^A>'¥^jj^j^(f) (4.11) 

k — 

^[K-i+u/iW - [log)"' + ¥(fc + 1 ) + ¥(fc + n/2)] 1+„/2(^) j 

Notice that the singular part (which is a polynomial in inverse powers of Vo") 
and the non-analytical parts (proportional to ^/o- and log cr) are expressed in terms 
of the the values of the function bq{A) at the integer points q, which are uniquely 
locally computable from the recursion relation, whereas the regular analytical part 
contains the values of the function bq{A) at half-integer positive points q and the 
derivatives of the function bq(A) with respect to q at integer positive points q, 
which are not expressible in terms of the local information. These objects are 
global and cannot be expressed further in terms of the local heat kernel coeffi- 
cients. However, they can be computed from the eqs. (|3.51l) and (13.541) in terms 
of the value of the function b(q) at some fixed point qo (see [3 J). 

The regular part of the Green function has a well defined diagonal value and 
the functional trace. It reads in odd dimensions n: 

Tr G''H^) = (-l)(«+i)/2(4;r)-"/2^L^5„/2-i(^) 

r(n/2) 
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(-ir.V^(4.,-.gl=|5;:iA. (4.12, 



and in even dimensions n 



Tr G'-^g(^) = {-\fl^-^^^^iB'^i^_,{A) - [log t^+^>{nll) - C] 5„/2-i(i)] 



r(n/2) 
(_l)«/2-i(4^)-n/2j 



^_^^„/2-l-/t 

)t!r(n/2-yt) 



A:=0 



C-W2-^) + log|-ij 



X ^-4] {-Xfl^-'-'Y{k+\-nl2)AA. (4.13) 



This trace determines the regularized vacuum expectation values like {ip^^ in 
quantum field theory. 
Thus, we see that 

i) all the singularities of the Green function and the non-analytical parts thereof 
(proportional to yfa in odd dimensions and to log cr in even dimensions) are 
determined by the values of the function bq{A) at integer points q, which are 
determined, in turn, by the heat kernel coefficients a^; 

ii) there are no power singularities, i.e. G^^"8(^) = q, in lower dimensions n = 
1,2; 

iii) there is no logarithmic singularity (more generally, no logarithmic part at 
all) in odd dimensions; 

iv) the regular part depends on the values of the function bq(A) at half-integer 
points q and its derivative bq(A) at integer points q and is a global object 
that cannot be reduced to purely local information like the heat kernel coef- 
ficients ak. 

The logarithmic part of the Green function is very important. On the one hand 
it determines, as usual, the renormalization properties of the regular part of the 
Green function, i.e. the derivative fi(d/dfi)G^^^(A). In particular. 



H^TrG^'H^) = \ 

Ofi 



for odd n 



I nn/2) 



■Bn/2-iW for even n. 
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On the other hand, it is of crucial importance in studying the Huygens prin- 
ciple. Namely, the absence of the logarithmic part of the Green function is a 
necessary and sufficient condition for the validity of the Huygens principle for 
hyperbolic operators. The heat kernel coefficients coefficients and, therefore, the 
logarithmic part of the Green function are defined for the hyperbolic operators just 
by analytic continuation from the elliptic case. Thus, the condition of the validity 
of Huygens principle reads 



or, by using (13.471 ). 



oo k-l+n/2 



kljir{k-j + n/2)\2) 

By expanding this equation in covariant Taylor series using the methods of [|3l one 
can obtain an infinite set of local conditions for validity of the Huygens principle, 
see [HI. In particular, 

[bn/2-lWf'^=0, (4.17) 

[V^K/2-i(A)t''^ = 0, (4.18) 

[V(^V,)ZP„/2-l(i)]'^'^S + ^g^,y[bn/2Wf'^ =0. (4.19) 



5 Heat Kernel Coefficients 

As we have shown above the calculation of the efi'ective action and the Green 
function reduces to the calculation of the heat kernel. An important part of that 
calculation is the calculation of the coefficients of the asymptotic expansion of the 
heat kernel. They are determined by a recursion system which is obtained simply 
by restricting the complex variable q in the eq. (13.571 ) to positive integer values 
q = 1,2,.... 
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5.1 Non-recursive Solution of the Recursion Relations 

This problem was solved in [|2l[3l[T0l where a systematic technique for calculation 
of ak was developed. The formal solution of this recursion system is 

1 \-'-/ 1 r'- / 1 



"'^ri^j ^v^—^n i-i>+T^j (5.1) 

Now, the problem is to give a precise practical meaning to this formal operator 
solution. To do this one has, first of all, to define the inverse operator (1 +D/ky^. 
This can be done by constructing the complete set of eigenvectors of the oper- 
ator D. However, first we introduce some auxiliary notions from the theory of 
symmetric tensors. 

Let 5"j be the bundle of symmetric tensors of type (m,n). First of all, we define 
the exterior symmetric tensor product 



V: S^^xS^^S';::^ (5.2) 



of symmetric tensors by 



(A V Bt\-^"^' . = £^"+1-^"+'^ . (5.3) 

V jai...am+j (ai...a,„ a,„+[...a„+j) ^ 

Next, we define the inner product 

by 

(A ★ Bfa\,,^a„, = ^a\.7a,„Bj[...y„ ■ (5.5) 

Finally, let I(„) be the identity endomorphism on the space of symmetric tensors 
of type (n,0); it is a section of the bundle S", that is, 

Ult=^-€^- (5-6) 



We also define the exterior symmetric covariant derivative 

' m 
'n+l 



V^: 5,T^ 5,7^1 (5.7) 



by 

(V^A)^;;;tv,=V(«,Af-^';;^^3. (5.8) 

These definitions are naturally extended to End ("V)- valued symmetric tensors, i.e. 
to the sections of the bundle S"^ ^EndCV). 
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5.2 Covariant Taylor Basis 

Let us consider the space of smooth two-point functions in a small neighborhood 
of the diagonal jc = jc'; we will denote such functions by |/). Let us define a special 
set of such functions { |n)}^Q, labeled by a non-negative integer n, by 

10) = 1, (5.9) 

\n) = V' (5.10) 

n\ 

where are the geometric parameters (normal coordinates) defined by (13.331) . 
These functions are scalars at the point x and symmetric tensors of type (0, n) at 
the point x' . It is easy to show that these functions satisfy the equation 

D\n)=n\n), (5.11) 

and, hence, are the eigenfunctions of the operator D with positive integer eigen- 
values. 

Let {n\ denote the dual linear functional defined by 

<n|/) = (V^)'7 (5.12) 

so that 

{n\m) = 5mnl(n) , (5.13) 

Using this notation the covariant Taylor series for an analytic function |/) can be 
written in the form 

oo 

\f) = J]\n)i^{n\f), (5.14) 

n=0 

For smooth functions the Taylor series is only an asymptotic series, which 
does not necessarily converge. For analytic functions, however, the Taylor series 
converges in a sufficiently small neighborhood of the fixed point x' . Therefore, 
the functions \n) form a complete orthonormal basis in the subspace of analytic 
functions. This is a reflection of the fact that an analytic function that is orthogonal 
to all functions \n), that is, whose all symmetrized derivatives vanish at the point 
x', is, in fact, identically equal to zero in this neighborhood. Note, however, that 
the space of functions we are talking about is not a Hilbert space since there are 
many analytic functions |/) such that the norm (/|/) defined above diverges. If 
we restrict ourselves to polynomials of some order, then this problem does not 
appear, and, hence, the space of polynomials is a Hilbert space with the inner 
product defined above. 
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5.3 Matrix Algorithm 

The complete set of eigenfunctions \n) can be employed to present the action of 
the operator L on a function |/) in the form 

Uf)= Yj \fn)'^iMUn)'k{n\f), (5.15) 

m,n>0 

where {m\L\n) are the 'matrix elements' of the operator L that are just EndCV)- 
valued symmetric tensors, i.e. sections of the vector bundle (8>End('y). When 
acting on an analytic function this series is nothing but the Taylor series of the 
result and converges in a sufficiently small neighborhood of the point x'\ for a 
smooth functions it gives an asymptotic epxansion. 

Now it should be clear that the inverse of the operator (l + jD^ in can be 
defined by 

(1 °° k 

l/>=ElT^I">*<"l/>- (5.16) 
' n=Q 

Using such representations for the operators (l + jD^ ^ and L we obtain a covari- 
ant Taylor series for the coefficients ak 

oo 

ak = Y^\n) -k {n\ak) (5.17) 



n=0 



where 



{n\ak) = y 

^ ^k + n k-l+nk-i 2 + n2 l+ni 

ni,...,nk-i>0 

X<n|L|nfc_i) ★ {nk-i\L\nk-2) ★ • ■ ■ ★ {ni\L\0) , (5.18) 

where the summation is over all non-negative integers ni,.. . ,nyt-i- It is not dif- 
ficult to show that for a differential operator of second order the matrix elements 
{m\L\n) do not vanish only for n < m + 2. Therefore, the summation over here is 
limited from above by 

0<ni, ni<ni+i+2, {i = l,2,...,k-l), (5.19) 



where n^ = n. Thus, the sum (15.181) contains only a finite number of terms. 
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Thus, we have reduced the problem of computation of the heat kernel coeffi- 
cients ak to the evaluation of the matrix elements {m\L\n) of the operator L. For 
a differential operator L of second order, the matrix elements {m\L\n) vanish for 
n > m + 2. Therefore, the summation over in (|5.18l) is limited from above: 
ni > 0, and ni < +2, for i = 1,2, ...,k- 1, and, hence, the sum (15.181) always 
contains only a finite number of terms. 

The matrix elements (n\L\nij of a Laplace type operator have been computed 
in our papers [131 [JJ. They have the following general form 

{m\L\m + 2) = -/vl(„), (5.20) 
<m|L|m+l) = 0, (5.21) 

(m\ ( ^ \ ( ^ \ 

(5.22) 

where g* is the metric on the cotangent bundle, Z(„) is a section of the vector 
bundle 5„ ^EndCV), Y(„) is a section of the vector bundle 5^ (S)End('y) and 
is a section of the vector bundle 5*^ (a symmetric tensor of type (2,n)). Here it is 
also meant that the binomial coefficient (^^ is equal to zero if k <0 or n <k. 

We will not present here explicit formulas, (they have been computed explic- 
itly for arbitrary m, n in [3,. lOJ), but note that all these quantities are expressed 
polynomially in terms of three sorts of geometric data: 

i) symmetric tensors of type (2,n), i.e. sections of the bundle Sj^ obtained by 
symmetric derivatives 

K(^n) = (y^T'^Riem (5.23) 



of the symmetrized Riemann tensor Riem taken as a section of the bundle 

J 2' 

ii) sections 

'R(n) = (y^r-''R (5.24) 



of the vector bundle 5/, <8iEnd('y) obtained by symmetrized derivatives of 
the curval 
End(^), 



the curvature 7? of the connection taken as a section of the bundle S | 
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iii) End CV)- valued symmetric forms, i.e. sections of the vector bundle S„<S) 
EndCV), constructed from the symmetrized covariant derivatives 

e(n) = (v^re (5.25) 

of the endomorphism Q. 

From dimensional arguments it is obvious that the matrix elements {n\L\n) are 
expressed in terms of the Riemann curvature tensor, Riem, the bundle curvature, 
K, and the endomorphism Q\ the matrix elements (n + 1 \L\n) — in terms of the 
quantities VRiem, VR and VQ; the elements {n + 2\L\n) — in terms of the quanti- 
ties of the form VVRiem, Riem ■ Riem, etc. 

5.4 Diagramatic Technique 

In the computation of the heat kernel coefficients by means of the matrix algorithm 
a "diagrammatic" technique, i.e., a graphic method for enumerating the difl^erent 
terms of the sum (15.181 ). turns out to be very convenient and pictorial [|3llT0l. 

The matrix elements {m\L\n) are presented by some blocks with m lines com- 
ing in from the left and n lines going out to the right (Fig. 1), 




Fig. 1 

and the product of the matrix elements {m\L\k) ★ {k\L\n) — by two blocks con- 
nected by k intermediate lines (Fig. 2), 




Fig. 2 

that represents the contractions of the corresponding tensor indices (the inner 
product). 

To obtain the coefficient {n\ak) one should draw, first, all possible diagrams 
which have n lines incoming from the left and which are constructed from k blocks 



Ivan G. Avramidi: Effective Action in Quantum Gravity 



41 



connected in all possible ways by any number of intermediate lines. When doing 
this, one should keep in mind that the number of the lines, going out of any block, 
cannot be greater than the number of the lines, coming in, by more than two and 
by exactly one. Then one should sum up all diagrams with the weight determined 
for each diagram by the number of intermediate lines from the analytical formula 
(15.181 ). Drawing of such diagrams is of no difficulties. This helps to keep under 
control the whole variety of different terms. Therefore, the main problem is re- 
duced to the computation of some standard blocks, which can be computed once 
and for all. 

For example, the diagrams for the diagonal values of the HMDS-coefficients 
^r^^ = (0|ayt) have the form. 



a^^= O (5.26) 



af^ = O O + ^ OX) (5.27) 



af^ 



o o o +^ o OX) +\ or) o (5.28) 

2 1 2 1 2 1 

+4-2 cixx) +4-3 cmx) +4-5 ■ 

As an illustration let us compute the coefficients a^^^^ and a^^^. We have ^ 

O =(0|LIO)=Z(o) (5.29) 

= (0ILI2) = (5.30) 

^ =<2ILI0)=Z(2)ai (5.31) 

CXI3 = (0|L|2) ★ <2|L|0) = -g'^^Dab , (5.32) 



where 



Z(0) = e-7^Ii/, (5.33) 
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Zi2)ab = '7(a'^b)Q-\'Rc(a'R'b) + ^^(a^\c\'R%) (5.34) 

Here, as usual, the parenthesis denote the complete symmetrization over all in- 
dices included and the vertical lines indicate the indices excluded from the sym- 
metrization. Hence, we immediately get 

af^ = Q-\Rlv, (5.35) 

and, by taking the trace of Z(2) and using the identity VaV^T?"^ = 0, we obtain the 
well known result [3J 

diag _ lpf,,\ _1a/0 , l<i? , ¥.,/j_AP_ J_p .pa^' , J_p , nabcd\ 



«2 =\Q- l^^y] - + -^KabK"" + IV ^Y^Ai? - ^RabR"" + ^RabcdR' ^ 

(5.36) 

The technique described above is manifestly covariant and is applicable for 
any Riemannian (or pseudo-Riemannian) manifold M and for any vector bundle 
V. It is also valid for local analysis on noncompact manifolds and manifolds with 
boundary (on a finite distance from the boundary). This method gives not only 
the diagonal values of the heat kernel coefficients but also the diagonal values 
of all their derivatives, that is, it gives also the off-diagonal coefficients in form 
of a covariant Taylor series. Due to the use of symmetric forms and symmetric 
covariant derivatives the famous 'combinatorial explosion' in the complexity of 
the heat kernel coefficients is avoided. This technique is very algorithmic and 
well suited to automated computation. The developed method is very powerful; it 
enabled us to compute for the first time the diagonal value of the fourth HMDS- 
coefficient ^ It was used in [331 [3T]| to compute the coefficient a^^^. 

Lastly, this technique enables one not only to carry out explicit computations, but 
also to analyse the general structure of the heat kernel coefficients for all orders k. 



5.5 General Structure of Heat Kernel Coefficients 

Now we are going to investigate the general structure of the heat kernel coeffi- 
cients. We will follow mainly our papers [fT0l l9l[TT]|. 

Our analysis will be again purely local. Since locally one can always expand 
the metric, the connection and the endomorphism Q in the covariant Taylor series, 
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they are completely characterized by their Taylor coefficients, i.e. the covariant 
derivatives of the curvatures, more precisely by the objects and 2(„) 

introduced above. We introduce the following notation for all of them 

5l(„) = {i^(„+2),!^(„+i),e(„)}, (n = 0,1,2,...), (5.37) 

and call these objects covariant jets; n will be called the order of a jet It is 
worth noting that the jets are defined by symmetrized covariant derivatives. This 
makes them well defined as the order of the derivatives becomes not important. It 
is only the number of derivatives that plays a role. 

The low-order coefficients Aq and Ai have been described above. As far as 
the higher order coefficients Ak, (k > 2), are concerned they are integrals of local 
invariants which are polynomial in the jets. One can classify all the terms in 
them according to the number of the jets and their order. The terms linear in 
the jets in higher order coefficients A^, (k > 2), are given by integrals of total 
derivatives, symbolically J^dvol tryA^^^%. They are calculated explicitly in ]^ 
\TOi\V\. Since the total derivative do not contribute to an integral over a complete 
compact manifold, it is clear that the linear terms vanish. Thus Ak, {k = 2,3,.. .), 
begin with the terms quadratic in the jets. These terms contain the jets of highest 
order (or the leading derivatives of the curvatures) and can be shown to be of the 
form J^dvol try^lA'^"^^. Then it follows a class of terms cubic in the jets etc. 
The last class of terms does not contain any covariant derivatives at all but only 
the powers of the curvatures. In other words, the higher order HMDS-coefficients 
have a general structure, which can be presented symbolically in the form 

Thus, for > 2 one can classify the terms in A2k according to the number of 
the jets and their order 

k 

Ak = Y,Ak,u), (5.38) 

./=2 

where Akxj) is the contribution of order j in the jets; they can be presented sym- 
bolically in the form 

Ak,(2) = Jjvoltrv^ 51(0)^(2^-4), (5.39) 

M 

^ 2k-6 

Ak,i3) = j dwo\trvJ]J]%iO)%ii)%(2k-6-i), (5.40) 
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,(5.41) 



M 

Akxk) = Jdwo\tvvj](^(0)f- (5-42) 



M 



More precisely, the functionals Akxj) transform under the rescaling of the jets 

%^k) ^ sa'^^k) (5.43) 

as follows 

A^^Oj^e-'a^Cfc-^Xo-)- (5-44) 

6 High Energy Approximation 

One can show that all quadratic terms can be reduced to five independent invari- 
ants, viz. ^B\M 

A,,(2) = ||^Jrfvoltr,|4i)eA^-2e + 2/f7?^^V,A'^-3v,7?«, 

M 

+/f ^ QA'^-^R + 4%«fcA^-2i?«* + /f >i?A^-27?|, (6.1) 

where /^'^ are some numerical coefficients. These numerical coefficients can be 
computed by the technique developed in the previous section. From the formula 
(15.181) we have for the diagonal coefficients a^^^^ up to cubic terms in the jets 

af'^ = {0\ak) = ^^^^{0;k-mO) 

V k ) 

k-l lik-i-l) p^~l) 

!=1 n,=0 [ k j\ i I 



where 



+0{%\ (6.2) 
<n; k\L\m) = (V V) * {n\L\m) (6.3) 
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and 0{%^) denote terms of third order in the jets. 

By computing the matrix elements in the second order in the jets and integrat- 
ing over M one obtains 111 |3l 

= 1' (6-4) 



1 



k-l 



/r = ^ ' (6-7) 

2(4A:2-1)' 

fi' = n ■ (6.8) 

4(4A:2-1) ^ ^ 

One should note that the same results were obtained by a completely different 
method in IfTTH . 

Let us consider the situation when the curvatures are small but rapidly varying 
(high energy approximation in quantum field theory), i.e. the derivatives of the 
curvatures are more important than the powers of them. This corresponds to an 
asymptotic expansion in the deformation parameter e as e ^ 0. Then the leading 
derivative terms in the heat kernel are the largest ones. Thus the heat trace has the 
form ^ 

0(0 ~ {Antr'l^ |ao - ?Ai + y^2(0| + 0{%\ (6.9) 

where H2{t) is some complicated nonlocal functional that has the following asymp- 
totic expansion as ? ^ 

H2(t)-22]-^AkX2). (6.10) 

k=2 

Using the results for Akxi) one can easily construct such a functional H2 just by a 
formal summation of the leading derivatives 

H2{t) = Jdvo\trvi^Qy^^\-tA)Q + 2'R\Va^y^^\-tA)Vb'R''' 

M 

-2Qy^^\-tA)R + Raby^^\-tA)R"^+Ry^^\-tA)R\, (6.1 1) 



1 
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where 7^'^(z) are entire functions defined by [[BO 
where 



= 


1, 


(6.13) 
(6.14) 


= 




(6.15) 


= 


6^ ' 


(6.16) 


f^'\o = 


48 ^ ' 


(6.17) 



Therefore, H2{t) can be regarded as generating functional for quadratic terms 
^yt,(2) (leading derivative terms) in all coefficients A^. It also plays a very impor- 
tant role in investigating the nonlocal structure of the effective action in quantum 
field theory in high-energy approximation [[BEl. 



7 Low Energy Approximation 

Let us consider now the opposite case, when the curvatures are strong but slowly 
varying (low-energy approximation in quantum field theory), i.e. the powers of 
the curvatures are more important than the derivatives of them. This corresponds 
to the asynptotic expansion in the deformation parameter or as a ^ 0. The main 
terms in this approximation are the terms without any covariant derivatives of the 
curvatures, i.e. the lowest order jets. We will consider mostly the zeroth order of 
this approximation which corresponds simply to covariantly constant background 
curvatures 

VRiem = 0, V:^ = 0, VQ = 0. (7.1) 
The asymptotic expansion of the heat trace 

0(0 ~ {Antr'i^Yj—r^k'^^> (^-^^ 
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determines then all the terms without covariant derivatives (highest order terms 
in the jets), (fc), in all heat kernel coefficients A^. These terms do not contain 
any covariant derivatives and are just polynomials in the curvatures and the en- 
domorphism Q. Thus the heat trace is a generating functional for all heat kernel 
coefficients for a covariantly constant background, in particular, for all symmetric 
spaces. Thus the problem is to calculate the heat trace for covariantly constant 
background. 



7.1 Algebraic Approach 

There exist a very elegant indirect way to construct the heat kernel without solving 
the heat equation but using only the commutation relations of some covariant first 
order differential operators. Below we wil follow our papers flU [51 [6l 121 [I2l [T3ll . 
The main idea is to employ a generalization of the usual Fourier transform to the 
case of operators; it consists in the following. We are going to use the following 
representation of the heat trace 

0(0 = ^ d\o\ ixv{Qx^{-tL)6{x,x')f'''^ . (7.3) 

M 

Let us consider for a moment a trivial case, where the curvatures vanish but 
the potential term does not: 

Riem = 0, 5^ = 0, VQ = 0. (7.4) 

In this case the operators of covariant derivatives obviously commute and form 
together with the potential term an Abelian Lie algebra 

[V^,Vv]=0, [V^,e]=0. (7.5) 

It is easy to show that the heat semigroup can be presented in the form 

exp(-?L) = (4;rO""^^exp(-/'0 J dk g^/^expl^-^^{k,gk) + k-V^, (7.6) 

where {k,gk) = k^gjjyk^ and k-V = k^Vj_i. Here, of course, it is assumed that the 
covariant derivatives also commute with the metric 



[V,^]=0. 



(7.7) 
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Acting with this operator on the Dirac distribution and using the obvious relation 



Of course, on curved manifolds the covariant differential operators V do not 
commute — their commutators are determined by the curvatures 21. The commuta- 
tors of covariant derivatives V with the curvatures % give the first derivatives of the 
curvatures, i.e. the jets ^(i), the commutators of covariant derivatives with 
give the second jets 21(2), etc. Thus the operators V together with the whole set of 
the jets J' form an infinite dimensional Lie algebra Q = {V, 2l(/); (/ = 1,2, . . .)}. 

Now, let us remember that the heat trace is a functional of the jets, with the 
jets being defined by symmetrized covariant derivatives. This makes the order of 
a jet well defined. For example, the structures involving commutators of covariant 
derivatives, (like [V^, Vh]R^/ ^, which involve 2-jets of the Riemann tensor on the 
left but, after using the Ricci identity, only 0-jets on the right) are not allowed. 
After symmetrizing over abed this jet vanishes. So, if we express the final answer 
for the heat kernel diagonal or for the heat kernel coeflicients in terms of the 
symmetrized jets, then there is a natural filtration with respect to the order of the 
jets involved. In other words, one can always say, what is the maximal order of 
symmetrized covariant derivative of the curvature involved in the result. This is 
especially true for the heat kernel coefficients Ak since they are polynomial in the 
jets. 

If we identify a small deformation parameter a with each derivative then a 
jet of order n is, actually, of order a". Thus, we get a perturbation theory in this 
small parameter. Since the derivatives are naturally identified with the momen- 
tum (or energy), the physicists call a situation when the derivatives are small the 
low-energy approximation. To evaluate the heat kernel in the low-energy approx- 
imation one can take into account only a finite number of low-order jets, i.e. the 
low-order covariant derivatives of the background fields, {2l(;); (i < N)} with some 
fixed A^, and neglect all higher order jets, i.e. the covariant derivatives of higher 
orders, i.e. put 2l(,) = for i > N. Then one can show that there exist a set of 
covariant differential operators that together with the low-order jets generate a 
finite-dimensional Lie algebra Qj^j = {V,2l(;);(/ = 1,2, . . . ,A^)}. One should stress 



[exp(A:- V)5(jc,jc')] = S(k), 



(7.8) 



one can integrate easily over k to obtain the heat trace 




(7.9) 



M 
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here what problem one can solve this way. We try to answer the following con- 
crete question: how do the heat kernel coefficients look if we through away all the 
(symmetrized) jets of order higher than A'^? 

Thus one can try to generalize the above idea in such a way that (|7.6I) would 
be the zeroth approximation in the commutators of the covariant derivatives, i.e. 
in the curvatures. Roughly speaking, we would like to find a representation of the 
heat semigroup in the form 

exp(-rL) = (Ant)-^'^ J dkO{t,k)exp^—^^{k,'i'{t)k) + T{k)j , (7.10) 

where {k,^{t)k) = k^"^ ABit)k'^ , T{k) = k^TA, (A = 1,2,...,D), are some first 
order differential operators and the functions ^{t) and 0(?,fc) are expressed in 
terms of commutators of these operators, i.e., in terms of the curvatures; that is, 
these functions are analytic functions of t. In general, the operators Ta do not 
form a closed finite dimensional algebra because at each step, by taking more 
commutators, there appear more and more derivatives of the curvatures. It is 
the low-energy reduction Q Qi^, i.e. the restriction to the low-order jets, that 
actually closes the algebra Q of the operators Ta and the background jets, i.e. 
makes it finite dimensional. 

Using this representation one can, as above, act with exp \T{ky\ on the Dirac 
distribution to get the heat kernel. The main point of this idea is that it is much 
easier to calculate the action of the exponential of the first order operator T{k) on 
the Dirac distribution than that of the exponential of the second order operator L. 



7.2 Covariantly Constant Background in Flat Space 

Let us consider now the more complicated case of nontrivial covariantly constant 
curvature of the connection on the vector bundle V in flat space: 

Riem = 0, Vn = 0, VQ = 0. (7.11) 

Using the condition of covariant constancy of the curvatures one can show that 
in this case the covariant derivatives form a nilpotent Lie algebra |i4j| 

[V^,Vv]=!??^v, (7.12) 
[V^,5^.^] = [V^,e]=0, (7.13) 
VR^v,Kp\ = VR^,,Q\=0. (7.14) 
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For this algebra one can prove a theorem expressing the heat semigroup oper- 
ator in terms of an average over the corresponding Lie group HI 



exp(-rL) = {Ant)-"'^ ex^{-tQ) 



det 



TM 



tn 



sinh(Z^) 



1/2 



X 



^ dk ^i/2exp|-i- {k,gt'RcotHt'R)k) + k-V 



(7.15) 
(7.16) 



where ^ • V = fc^V^. Here functions of the curvatures are understood as functions 
of sections of the bundle End(rM)(8)End('y), and the determinant detj-M is taken 
with respect to the tangent space indices; the fiber indices of the bundle "V being 
intact. 

It is not difficult to show that in this case we also have 



[exp(k-V)6{x,x)f''^ = 6{k). 



(7.17) 



Subsequently, the integral over k becomes trivial and we obtain immediately the 
trace of the heat kernel [|4l 



0(0 



= (4nt)-"'^ J 



dvol trvexp(-?0 



det 



TM 



tn 



sinh(?7?) 



1/2 



(7.18) 



M 



Expanding it in a power series in t one can find all covariantly constant terms in 
all heat kernel coefficients Ak. 

As we have seen the contribution of the curvature "^^y is not as trivial as that of 
the potential term. However, the algebraic approach does work in this case too. It 
is a good example how one can get the heat kernel without solving any differential 
equations but using only the algebraic properties of the covariant derivatives. 



7.2.1 Quadratic Potential in Flat Space 

In fact, in flat space it is possible to do a bit more, i.e. to calculate the contribution 
of the first and the second derivatives of the potential term Q [6J. That is, we 
consider the case when the derivatives of the endomorphism Q vanish only starting 
from the third order, i.e. 



Riem = 0, Vn = 0, VVVQ = 0. 



(7.19) 
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Besides we assume the background to be Abelian, i.e. all the nonvanishing back- 
ground quantities, %ap, Q, Q;^i = ^jjQ and Q-y^j = V^VyQ, commute with each 
other. Thus we have again a nilpotent Lie algebra 

[V^,Vv]=!??^v, (7.20) 
[V^,e] = e;^, (7.21) 

[V^,e;v] = e;v^, (7.22) 

all other commutators being zero. 

Now, let us represent the endomorphism Q in the form 

Q = Qo-a'''NiNk, (7.23) 

where (z = I,..., q;q <n), a'^ is some constant symmetric nondegenerate qxq ma- 
trix, Qo is a covariantly constant endomorphism and Ni are some endomorphisms 
with vanishing second covariant derivative: 

VQo = 0, VVNi = 0. (7.24) 

Next, let us introduce the operators Xa = (V^,A^,), (A = l,...,n + q) and the matrix 

irAB) = ('^// ^i:'], (7.25) 



-Nk-y 

withA^/;^ = V^A^,-. 

The operator L can now be written in the form 

L = -G^''XaXb + Qo, (7.26) 

where 



(C) = (o ^k)- (7.27) 

and the commutation relations (|7.22l) take a more compact form 

[Xa,Xb]=Tab, (7.28) 

all other commutators being zero. 

This algebra is again a nilpotent Lie algebra. Thus one can apply the previous 
results in this case too to get [6J 

1/2 

exp(-?L) = (47rO~^"^'^^/^exp(-?2o) 



det 



tT 



sinh(r:r) / 



X ^ dkG^I^Qx^{^j^{k,GtTcoi\\{t^)k) + X{k)^, (7.29) 
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where G = detGAs and X{k) = R^Xa- 

Thus we have expressed the heat semigroup operator in terms of the operator 
exp {X{k)\ . The integration over k is Gaussian except for the noncommutative part. 
Splitting the integration variables {k^) = {q^,a)') and using the Campbell-Hausdorf 
formula we obtain ^ 

[exp [X(k)] 6{x, x')f'''^ = exp [N{oj)] 6(q), (7.30) 

where N{oj) = a/Ni. Further, after taking off the trivial integration over q and a 
Gaussian integral over co, we obtain the heat trace 161 . 

To describe the result let us introduce a matrix determined by second deriva- 
tives of the potential term as follows 

i' = (p^), /'^ = ^v^v,e, (7.31) 

and the matrices C{t) = (C^v(O), Kit) = (K^'yit)) S (t) = (5^v(0) and Eit) = (E^'yit)) 
by 

C(t) = , (7.32) 
c 

m = ^^F(z)^sinh(^) , (7.33) 

c 

5(0 = ^^TO^sinh(^), (7.34) 

c 

E{t) = j)^F{z)t%m\ii^-^, (7.35) 



where 



F{z) = {l-7^-z^P)-^ (7.36) 



and the integral is taken along a sufficiently small closed contour C that encircles 
the origin counter-clockwise, so that F{z) is analytic inside this contour. 
Then the heat trace has the form 



0(0 = (Ant)-"'^ J dwol try [0(0]"^/^ exp 



M 



-tQ + ^t\VQ,l'{t)VQ) 



(7.37) 
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where {VQ,^{t)VQ) 



det xMK{t) det tm 

[l+t^Cit)P] 
xdetTM{l+t^[E(t)-S{t)K-\t)S(t)]P] , 

(^^y(t)) = [l+t^C(t)P]~'c{t). 



(7.38) 



T(0 



(7.39) 



The formula (|7.37l) exhibits the general structure of the heat trace. One sees 
immediately how the endomorphism Q and its first derivatives VQ enter the re- 
sult. The nontrivial information is contained only in a scalar, 0(?), and a ten- 
sor, ^^v(0- These objects are constructed purely from the curvature "^^y and the 
second derivatives of the endomorphism Q, WQ. Thus, the heat kernel coeffi- 
cients Ak are constructed from three difi'erent types of scalar (connected) blocks, 
Q, 0(„)(7?, VV(2) and V^2*Pj';^(^, VVe)Vy<2. They are listed explicitly up to Ag 
inflH. 

7.3 Homogeneous Bundles over Symmetric Spaces 

The exposition in this section closely follows our papers ^ U\ [IS Il3l. Our goal 
is to compute the heat kernel of the Laplace type operator L = -A -I- 2 in the zero- 
order of the low-energy approximation. The difference with the previous sections 
is that now we are going to do it on most general covariantly constant background, 
that is, bundles with parallel curvature (that are called homogeneous bundles) on 
Riemannian manifolds with parallel curvature (that are called symmetric spaces). 

It is well known that heat invariants are determined essentially by local ge- 
ometry. They are polynomial invariants in the curvature with universal constants 
that do not depend on the global properties of the manifold [27|. It is this univer- 
sal structure that we are interested in this paper. Our goal is to compute the heat 
kernel asymptotics of the Laplacian acting on homogeneous vector bundles over 
symmetric spaces. 

In this section we will further assume that M is a locally symmetric space with 
a Riemannian metric with the parallel curvature 



which means, in particular, that the curvature satisfies the integrability constraints 



(7.40) 




bed - R-^^ ebR'^ acd + ec-R^ dab - edR'^ cab = 0. 



(7.41) 
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In the following we will also consider homogeneous vector bundles with parallel 
bundle curvature 

V;,^a/j = 0, (7.42) 
which means that the curvature satisfies the integrability constraints 

VRcd,'Rab\ -R^acd'Rfb -R^hcdKf = 0. (7.43) 

Finally, we consider a parallel section Q of the endomorphism bundle End('y), 
that is, 

VpQ = 0, (7.44) 

which means that 

[Kd,Q]=0. (7.45) 

We will use normal coordinates defined above. Note that for symmetric spaces 
normal coordinates cover the whole manifold except for a set of measure zero 
where they become singular [|T9l . This set is precisely the set of points conjugate 
to the fixed point x' (where A~^(x,x') = 0) and of points that can be connected to 
the point x' by multiple geodesies. In any case, this set is a set of measure zero 
and, as we will show below, it can be dealt with by some regularization technique. 
Thus, we will use the normal coordinates defined above for the whole manifold. 

7.3.1 Curvature Group of a Symmetric Space 

We assumed that the manifold M is locally symmetric. Since we also assume 
that it is simply connected and complete, it is a globally symmetric space (or sim- 
ply symmetric space). A symmetric space is said to be compact, non-compact 
or Euclidean if all sectional curvatures are positive, negative or zero. A generic 
symmetric space has the structure M = Mq x Ms , where Mq = R"" and Ms is a 
semi-simple symmetric space; it is a product of a compact symmetric space M+ 
and a non-sompact symmetric space A/_, Ms = M+ x M_ . Of course, the dimen- 
sions must satisfy the relation no + ns = n, where = dimM^. 

Let A2 be the vector space of 2-forms on M at a fixed point x' . It has the 
dimension dim A2 = n(n - l)/2, and the inner product in A2 is defined by 

{X,Y) = ^X,hY"'^. (7.46) 
The Riemann curvature tensor naturally defines the curvature operator 



Riem : A2 ^ A2 



(7.47) 
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by 

{KitmX)ab = Uiab"^X,d. (7.48) 

This operator is symmetric and has real eigenvalues which determine the principal 
sectional curvatures. Now, let Ker(Riem) and Im(Riem) be the kernel and the 
range of this operator and 

n{n - 1) 

p = dimlm(Riem) = — dimKer(Riem) . (7.49) 

Further, let Ai, (z = 1, ... ,p), be the non-zero eigenvalues, and E'at be the corre- 
sponding orthonormal eigen-two-forms. Then the components of the curvature 
tensor can be presented in the form [|7]| 

Rabcd=/3ikE'abE\d, (7.50) 

where fiik is a symmetric, in fact, diagonal, nondegenerate pxp matrix. Of course, 
the zero eigenvalues of the curvature operator correspond to the flat subspace Mq, 
the positive ones correspond to the compact submanifold M+ and the negative 
ones to the non-compact submanifold M_. Therefore, Im(Riem) = T^Ms- 

In the following the Latin indices from the middle of the alphabet will be used 
to denote tensors in Im(Riem); they should not be confused with the Latin indices 
from the beginning of the alphabet which denote tensors in M. They will be raised 
and lowered with the matrix /3ik and its inverse yS'^. 

Next, we define the traceless nxn matrices D,- = {D"ib), where 

D"ib = -fi,kE\-bS'". (7.51) 

The matrices Dj are known to be the generators of the holonomy algebra, 'K, i.e. 
the Lie algebra of the restricted holonomy group, H, 

[Di,Dk]=FjikDj, (7.52) 

where pJik are the structure constants of the holonomy group. The structure con- 
stants of the holonomy group define the pXp matrices Fj, by (F,)4 = F^ik, which 
generate the adjoint representation of the holonomy algebra, 

[Fi,Fk]=FJ,kFj. (7.53) 

For symmetric spaces the introduced quantities satisfy additional algebraic 
constraints. The most important consequence of the eq. (17.411) is the equation |I71 

E'ac-D^kb - E'bcD^ka = F'kjE-'ab ■ (7.54) 
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Now, we introduce a new type of indices, the capital Latin indices, A,B,C,..., 
which split according to A = (a,/) and run from 1 to N = p + n. We define new 
quantities C^bc by 

C'ab = E'ab, C"ib = -Chi = D'^ih, C'ki = F'ki, (7.55) 

all other components being zero. Let us also introduce rectangular pxn matrices 
Ta by {Tayc = E^ac and the n x p matrices Ta by (fa)^,- = -D^'ia. Then we can 
define NxN matrices Ca = (Ca, Ci) 

so that (Ca)^c = C^AC- 

Then one can prove the following [j7| : 

Theorem 1 The matrices Ca generate the adjoint representation of a Lie algebra 
Q with the structure constants C^bc, that is, 

[Ca,Cb] = C^abCc, (7.57) 

For the lack of a better name we call the algebra ^ the curvature algebra. As 
it will be clear from the next section it is a subalgebra of the total isometry algebra 
of the symmetric space. It should be clear that the holonomy algebra "K is the 
subalgebra of the curvature algebra The curvature algebra is compact; it is a 
direct sum of two ideals, Q = Qq ®Qs, an Abelian center Qq of dimension nQ and 
a semi-simple algebra Qs of dimension p + n^. 

Next, we define a symmetric nondegenerate NxN matrix 

<''-) = (*o' ftj- ('-^^ 

This matrix and its inverse y'^^ will be used to lower and to raise the capital Latin 
indices. 



7.3.2 Killing Vectors Fields 

We will use extensively the isometrics of the symmetric space M. We follow 
the approach developed in [|7l |3l [TOl [131. The generators of isometrics are the 
Killing vector fields ^. The set of all Killing vector fields forms a representation 
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of the isometry algebra, the Lie algebra of the isometry group of the manifold 
M. We define two subspaces of the isometry algebra. One subspace is formed by 
Killing vectors (called translations) satisfying the initial conditions Vp^^| = , 
and another subspace is formed by Killing vectors (called rotations) satisfying the 
initial conditions £A , = . 

^ \x=x' 

One can easily show that a basis of translations can be chosen as 

Pa = [<Kcoi<Kf (7.59) 

where K = {K"h) is a matrix defined by 

K\=R\.bd/y'^. (7.60) 
We can also show that the vector fields 

U = -D',ay"-^, (7.61) 

define p linearly independent rotations. By adding the trivial Killing vectors for 
flat subspaces we find that the number of independent rotations is p + n^ns + 
no(^o - l)/2 . We introduce the following notation {^a) = {Pa,U)- 

By using the explicit form of the Killing vector fields obtained above Id one 
can prove the following theorem. 

Theorem 2 The Killing vector fields form form a representation of the curva- 
ture algebra Q 

[U,^b\ = C^ab^c. a. 62) 

Notice that they do not generate the complete isometry algebra of the symmetric 
space M. The curvature algebra Q introduced in the previous section is a subalge- 
bra of the total isometry algebra. It is clear that the Killing vector fields form 
a representation of the holonomy algebra "K, which is the isotropy algebra of the 
semi-simple submanifold Ms, and a subalgebra of the total isotropy algebra of the 
symmetric space M. 

7.3.3 Homogeneous Vector Bundles 

Let h'^b be the projection to the subspace TxMs of the tangent space and 



IM 

q b = o h-h h 



(7.63) 
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be the projection tensor to the flat subspace R"" Since the curvature exists only in 
the semi-simple submanifold Ms, the components of the curvature tensor Rabcd, 
as well as the tensors E'ab, are non-zero only in the semi-simple subspace T^M^. 
Then 

Rabcdq^e = Rabq% = E\bq% = D^ibqK = D",bqa" = . (7.64) 

Equation (|7.43l) imposes strong constraints on the curvature of the homoge- 
neous bundle 'W. We define 

Sab = n^q'aq^ , Gab = nl^h%h\ , (7.65) 

so that 

K^=&ab + Sab. (7.66) 

Then, from eq. (17.431 ) we obtain 

[!Bab,Sad] = [Sab,&cd]=0, (7.67) 

and 

[£'cd,£'ab] -R-^acd&fb -R-^bcd&af = . (7.68) 

This means that Sab takes values in an Abelian ideal of the gauge algebra ^ym 
and Sab takes values in the holonomy algebra. More precisely, eq. (17.681) is only 
possible if the holonomy algebra "TY is an ideal of the gauge algebra ^ym- Thus, 
the gauge group Gym must have a subgroup ZxH, where Z is an Abelian group 
and H is the holonomy group. 

Let Xab be the generators of the orthogonal algebra SO{n) is some representa- 
tion X. Then the matrices T, = -jD"ihX^\, are the generators of the gauge algebra 
^YM realizing a representation T of the holonomy algebra "K. Next, we can show 
that the curvature of the homogeneous bundle is given by 

^If = -E'abT:+!Bab = \R'''abXcd + Sab. (7.69) 

Now, we consider the representation S of the orthogonal algebra SO{n) defin- 
ing the spin-tensor bundle T and define the matrices 

Gab = ^ab ®lx + h.® Xab ■ (7.70) 

Obviously, these matrices are the generators of the orthogonal algebra SO{n) in 
the product representation E(8)Z. Next, the matrices 7,- = -jD^HjH^a form a rep- 
resentation Y of the holonomy algebra and the matrices 

% = -^D'',bG^a (7.71) 
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are the generators of the holonomy algebra in the product representation = Y®T. 

Then the total curvature, that is, the commutator of covariant derivatives, (13.61) 
of a twisted spin-tensor bundle "V is 

Kb = -E\b%^Bab = )f\bGcd^^ab. (7.72) 

7.3.4 Twisted Lie Derivatives 

Let ^ be a section of a twisted homogeneous spin-tensor bundle T . Let be the 
basis of Killing vector fields. Then the covariant (or generalized, or twisted) Lie 
derivative of along is defined by 

La^ = I^A^V^ + ^A^-bC'^ay. (7.73) 

One can prove the theorem [[T2l[T3]| . 
Theorem 3 The operators Xa satisfy the commutation relations 

[£a,£b] = C'^ab£c + !Bab, (7.74) 

where 

The operators Xa form an algebra that is a direct sum of a nilpotent ideal and 
a semisimple algebra. For the lack of a better name we call this algebra gauged 
curvature algebra and denote it by ^gauge- 

Now, let us define the operator 

£2 = /^x^Xg (7.76) 
and the Casimir operator of the holonomy group 

^2^ (7.77) 

Then one can prove that [fT3l 

Theorem 4 The Laplacian A acting on sections of a twisted spin-tensor bundle 
•y over a symmetric space has the form 



(7.78) 
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7.3.5 Geometry of the Curvature Group 

Let Ggauge be the gauged curvature group and H be its holonomy subgroup. Both 
these groups have compact algebras. However, while the holonomy group is al- 
ways compact, the curvature group is, in general, a product of a nilpotent group. 
Go, and a semi-simple group, G.v, Ggauge = Go x Gs ■ The semi-simple group Gs is 
a product Gs = G+ x G_ of a compact G+ and a non-compact G- subgroups. 

Let be the basis Killing vectors, be the canonical coordinates on the 
curvature group G and ^{k) = k^^A- The canonical coordinates are exactly the 
normal coordinates on the group defined above. Let Ca be the generators of the 
curvature group in adjoint representation and C{k) = k^CA- 

Let X = (Xa^) be the matrix defined by 

X= ^ , (7.79) 

l-exp[-C(yt)] 

and Xa be the vector fields on the group G defined by 

^A=^A^Ji7- (7.80) 

Then one can show that [|T3l the vector fields Xa form a representation of the 
curvature algebra ^ 

[Xa,Xb] = C^abXc. (7.81) 

The vector fields Xa are nothing but the right-invariant vector fields. 

Since we will actually be working with the gauged curvature group, we intro- 
duce now the operators (covariant right-invariant vector fields) J a by 

JA=XA-^BABk\ (7.82) 
Then we show [fT3ll that the operators J a form the following algebra 

Ua,Jb\ = C^abJc+^ab. (7.83) 
Thus, the operators Ja form a representation of the gauged curvature algebra 

^gauge- 

Now, let Xa be the Lie derivatives and £,{k) = k^£,A- Then we find |fT3l 

Ja cxp[£{k)] = cxp[£{k)]£A ■ (7.84) 

Notice that J a are first order diff"erential operators with respect to k^, whereas Xa 
are first-order partial diff'erential operators with respect to the coordinates x acting 
on sections of the bundle 'V. 
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7.3.6 Heat Kernel on the Curvature Group 

Now, let us define the operator 

72 = j^^JaJb (7.85) 
and the invariant (scalar curvature of the curvature group) 

^g = -^/''C^adCV. (7.86) 

Then by using the properties of the right-invariant vector fields J a one can find 
the heat kernel of the operator on the curvature group Q [fT3l . 

Theorem 5 Let 0(?; k) be a function on the curvature group defined in canonical 
coordinates by 



^{t;k) = i4ntr^^^ 



. ^ /sinh[rS]\l-i/2r /sinh[C(fc)/2]\ 



-1/2 



xexp|-^ {k,ytScothitS)k) + ^Rct^ , (7.87) 

where {u,yv) = jabu^v^ is the inner product on the algebra Q. Then 0(?;fc) satis- 
fies the heat equation 

dt^ = J^^, (7.88) 



and the initial condition 



where y = detyAS. 



^{0\k) = y-^l^6{k), (7.89) 



In the following we will complexify the gauged curvature group in the follow- 
ing sense. We extend the canonical coordinates {k^) = (p",a>') to the whole com- 
plex Euclidean space C^. Then all group-theoretic functions introduced above 
become analytic functions of k^ possibly with some poles on the real section 
for compact groups. In fact, we replace the actual real slice R.'^ of with an 
//-dimensional subspace R^^^ in obtained by rotating the real section coun- 
terclockwise in by n/4-. That is, we replace each coordinate k^ by e'^^^k^. In 
the complex domain the group becomes non-compact. We call this procedure the 
decompactification. If the group is compact, or has a compact subgroup, then this 
plane will cover the original group infinitely many times. 
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Since the metric (jab) = diag(5a/^,/3,j) is not necessarily positive definite, (ac- 
tually, only the metric of the holonomy group /3,j is non-definite) we analytically 
continue the function 0(?; k) in the complex plane of t with a cut along the negative 
imaginary axis so that -njl < arg t < 3n/2. Thus, the function 0(?;^) defines an 
analytic function of t and k"^. For the purpose of the following exposition we shall 
consider t to be real negative, t <0. This is needed in order to make all integrals 
convergent and well defined and to be able to do the analytical continuation. 

As we will show below, the singularities occur only in the holonomy group. 
This means that there is no need to complexify the coordinates p". Thus, in the 
following we assume the coordinates p" to be real and the coordinates co' to be 
complex, more precisely, to take values in the p-dimensional subspace R.feg of 
obtained by rotating RP counterclockwise by n/A in C^^ That is, we have R^g = 

This procedure (that we call a regularization) with the nonstandard contour of 
integration is necessary for the convergence of the integrals below since we are 
treating both the compact and the non-compact symmetric spaces simultaneously. 
Remember, that, in general, the nondegenerate diagonal matrix /3ij is not positive 
definite. The space Rfgg is chosen in such a way to make the Gaussian exponent 
purely imaginary. Then the indefiniteness of the matrix /3 does not cause any prob- 
lems. Moreover, the integrand does not have any singularities on these contours. 
The convergence of the integral is guaranteed by the exponential growth of the 
sine for imaginary argument. These integrals can be computed then in the follow- 
ing way. The coordinates ojj corresponding to the compact directions are rotated 
further by another n/4 to imaginary axis and the coordinates co-^ corresponding to 
the non-compact directions are rotated back to the real axis. Then, for ? < all the 
integrals below are well defined and convergent and define an analytic function of 
nn a complex plane with a cut along the negative imaginary axis. 

7.3.7 Heat Trace 

Now, by using the heat kernel (17.871) of the operator on the curvature group 
obtained above, the relation (17.781) of the Laplacian and the operator and the 
property (17.841) one can find the following integral representation of the heat semi- 
group of the Laplace-type operator [ISJ. 

Theorem 6 Let L = -A + Q be the Laplace type operator acting on sections of a 
homogeneous twisted spin-tensor vector bundle over a symmetric space. Then the 
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heat semigroup exp(-?L) can be represented inform of an integral 

xexpj-^ {k,ytScoih{tS)k)^QX^{L{k)] . (7.90) 

The heat trace can be obtained by acting by the heat semigroup exp(-?L) on the 
delta-function, To be able to use this integral representation we need to compute 
the action of the isometrics exp[X(fc)] on the delta-function. 

Let a/ be the canonical coordinates on the holonomy group H and {k^) = 
(p",co') be the natural splitting of the canonical coordinates on the curvature group 
G. Then we can prove that [|T3]| 



[cxpimmx,/)] 



diag _ 



det 



/sinh[D(a;)/2]y ^ 



oxpmiomp), (7.91) 



where D(co) = co'Di and 'RiaS) = (0%. 

We implicitly assumed that there are no closed geodesies and that the equa- 
tion of closed orbits of isometrics has a unique solution. On compact symmetric 
spaces this is not true: there are infinitely many closed geodesies and infinitely 
many closed orbits of isometrics. However, these global solutions, which reflect 
the global topological structure of the manifold, will not afl^ect our local analysis. 
In particular, they do not aff'ect the asymptotics of the heat kernel. That is why, 
we have neglected them here. This is reflected in the fact that the Jacobian in 
(17.911 ) can become singular when the coordinates of the holonomy group a;' vary 
from -oo to oo. Note that the exact results for compact symmetric spaces can be 
obtained by an analytic continuation from the dual noncompact case when such 
closed geodesies are absent [1X91 . That is why we proposed above to complexify 
our holonomy group. If the coordinates o)' are complex taking values in the sub- 
space Rfeg defined above, then the equation of closed orbits should have a unique 
trivial solution and the Jacobian is an analytic function. It is worth stressing once 
again that the canonical coordinates cover the whole group except for a set of 
measure zero. Also a compact subgroup is covered infinitely many times. 
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Now by using the above lemmas and the theorem we can compute the heat 
trace. We define the invariant (scalar curvature of the holonomy group) 



Theorem 7 The heat trace of the operator L has the form 

/sinh(/'S)\T^^^ 



0(0 = (4, 



X 



X 



■nt)-"l^ J 

I 



dvol tri 



M 
dco 



detT-M 



\ tS 



1 



exp 



(7.92) 



det 



{4nt)P/^ 

( 



fi^'^ ^^P I - ^ (^^P^) \ cosh [:^(a;)] 



/ sinh[F(fa;)/2] ^^^^^ 
Fico)/2 



det 



TM 



sinh[D(a;)/2] 



D(co)/2 



-1/2 



(7.93) 



where /3 = detjS/,-, {a),/3a)) = fiijio'ojj and F(a)) = oj'Fi. 



This equation can be used now to generate all heat kernel coefficients for 
any locally symmetric space simply by expanding it in a power series in t. By 
using the standard Gaussian averages one can obtain now all heat kernel coeffi- 
cients in terms of traces of various contractions of the matrices D^u, and F^jk with 
the matrix jS'^. All these quantities are curvature invariants and can be expressed 
directly in terms of the Riemann tensor. 



8 Low Energy Effective Action in Quantum General 
Relativity 

We can apply now the obtained results for the heat trace to compute the low- 
energy one-loop effective action in quantum general relativity given by (12.711 ). In 
the Euclidean formulation we have 

Td) = i(logDetL-21ogDetF) , (8.1) 

which, in the zeta regularization takes the form 



r(i) = -^fc(0)-2^;(0)), 



(8.2) 
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where ^i{s) and ^f{s) are the zeta functions of the graviton operator L and the 
ghost operator F. Now, let us define the total zeta function by 

^gr{s) = ^i{s)-1^f{s). (8.3) 

Then the effective action is 

r(i) = -^^^^(0). (8.4) 
Next, by using the definition of the zeta function we obtain 

^gr{s) = ^Adt t'-^e'^®GR{t) , (8.5) 


where 

®GR{t) = &i{t)-2@F{t), (8.6) 

is a renormalization parameter introduced to preserve dimensions and 0/^(0 
and 0f(O are the heat traces of the operators L and F. Here /I is a sufficiently 
large negative infrared cutoff parameter introduced to regularize any infrared di- 
vergences which are present if the operators L and F have negative modes. The 
parameter A should be set to zero at the end of the calculations. 

Now, notice that both operators L and F are of Laplace type, that is, -A + Q, 
acting on pure tensor bundles; so, there is no Yang-Mills group here, ^ab = ^ab = 
Bab = 0. The operator L acts on the bundle T(2) = T*M^ T*M of symmetric two- 
tensors and the operator F acts on sections of the tangent bundle T(i) = TM. The 
potentials, Q, for both operators are obviously read off from their definition 

(Q(i)fb = -R''b, (8.7) 

{Qi2)l/' = -2R^%'\-26\R''^d)+Rcd8''' + -^8cdR''' 

-T-^^gcdg^'R + 6\cS\){R - 2A) . (8.8) 
(n-2) 

The generators of the orthogonal group S 0(n) in the vector representation are 



{^il)ab) d = '^S^[agb}d- 



(8.9) 
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The generators of the orthogonal group S 0(n) in the symmetric 2-tensor repre- 
sentation are 

(^(2)ab) c/^ = -4S^'lagb]idS^\) . (8.10) 

The generators of the holonomy group are 



and 



%), = -2Av/(i), 

which, in component language, reads 



and 

The Casimir operators are 
and 



d) ■ 



(8.11) 
(8.12) 

(8.13) 

(8.14) 
(8.15) 
(8.16) 



By using the results for the heat traces described above we obtain the total heat 
trace 



®GRit) = (Ant)-"'^ J dwol exp 1(^7? + ii?// 1 ? 

M 



8 



(8.17) 



X 



^reg 



X det^ 



/sinh[F(cj)/2] 



Fico)/2 



1/2 



det 



TM 



sinh[DM/2] 
D(oj)/2 



-1/2 



(8.18) 



where 



'¥(t;oj) = exp[-?(i?-2A)]tr7-(2)exp(?y(2))cosh[2D(a;) V/(i)] 

-2trrMexp(?V(i))cosh[D(w)] , (8.19) 
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and the matrices V(i) and V(2) are defined by 

[Vii))\ = 2R\, (8.20) 

(V(2)) ;'=^S'-\cR'^d)-Rcdg''' ^-^gcdR''' + -\^gcdg'"R. (8.21) 

^ n-2 in -2) 

One can go further and compute the function *P(?; oS) by finding the eigenvalues 
of the endomorphisms V^i) and V(2)- However, we will not do it here and leave 
the answer in the form (|8.19l) . By using the obtained heat trace one can compute 
now the zeta function and then the effective action. We would like to stress two 
points here. First of all, quantum general relativity is a non-renormalizable theory. 
Therefeore, even if one gets a final result via the zeta-regularization one should 
not take it too seriously. Secondly, our results for the heat kernel and, hence, 
for the effective action are essentially non-perturbative. They contain an infinite 
series of Feynmann diagrams and cannot be obtained in any perturbation theory. 
One could try now to use this result for the analysis of the ground state in quantum 
gravity. But this is a rather ambitious program for the future. 
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